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Based on recent advances in first-principles theory, we develop a general model of the band 
offset at metal/ferroelectric interfaces. We show that, depending on the polarization of the film, 
a pathological regime might occur where the metallic carriers populate the energy bands of the 
insulator, making it metallic. As the most common approximations of density functional theory are 
affected by a systematic underestimation of the fundamental band gap of insulators, this scenario 
is likely to be an artifact of the simulation. We provide a number of rigorous criteria, together with 
extensive practical examples, to systematically identify this problematic situation in the calculated 
electronic and structural properties of ferroelectric systems. We discuss our findings in the context 
of earlier literature studies, where the issues described in this work have often been overlooked. We 
also discuss formal analogies to the physics of polarity compensation at LaAlOa /SrTiOa interfaces, 
and suggest promising avenues for future research. 

PACS numbers: 71.15.-m, 73.61.-r, 77.55.-s, 77.80-e 
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I. INTRODUCTION 

Advances in oxide thin film growth techniques over 
the last ten years have led to the fabrication of many 
novel oxide-based metal-insulator heterostructures with a 
dizzying range of functionalities. Not only are the current 
technological limits of information storage density and 
speed being pushed forward by the use of, e.g., nanoscale 
ferroelectric memorieSf^"— but entirely new concepts in 
device applications are also emerging, in which the elec- 
trical and the magnetic degrees of freedom are both 
present within the same active element and strongly cou- 
pled/'^ Examples of this trend include thin-film capaci- 
tors^ strongly correlated field-efltect devices^ and mag- 
netic/ferroelectric tunnel junctionSjAS.^— 

Density functional theory (DFT) methods, either 
within the local density (LDA) or generalized gradi- 
ent (GGA) approximation, have been an invaluable tool 
in achieving a fundamental understanding of this class 
of systems particularly with recent developments 

which allow the application of finite electric fields to pe- 
riodic solids or layered heterostructuresji^— However, 
since this domain of research is relatively new, it is im- 
portant to identify, in addition to the virtues, also the 
limitations of DFT that are specific to metal/ferroelectric 
interfaces, and that when overlooked might lead to erro- 
neous physical conclusions. 

For most practical applications, a capacitor must be 
insulating to DC current; transmission of electrons via 
non-zero conductivity and/or direct tunneling (leakage) 
is generally an undesirable source of heating and power 
consumption. At the quantum mechanical level, the in- 
sulating properties of a capacitor are guaranteed by the 
presence of a dielectric film with a finite band gap at the 
Fermi level, where propagation of the metallic conduction 



electrons is forbidden. In the language of semiconductor 
physics, we can alternatively say that both Schottky bar- 
rier heights (SBH), respectively and (pp for electrons 
and holes, need to be positive for the device to behave as 
a capacitor. (By convention we assume that, if the Fermi 
level of the metal lies in the gap of the insulator, both 
(/)„ and (/)p are positive.) 

If, on the contrary, either (pp or 0„ is negative, in- 
jection of holes or electrons into the dielectric becomes 
energetically favorable and the device behaves instead as 
an Ohmic contact. Most importantly, at such a junc- 
tion there is necessarily (at thermodynamic equilibrium) 
a spill-out of charge from the metal to the insulator, as 
the system re-equilibrates the chemical potential of the 
free carriers on either side. Such intrinsic space charge 
induces metallicity (by intrinsic doping) in the dielec- 
tric film, and overall profoundly alters the electronic and 
structural properties of the interface. 

While in principle the charge spillage might be a real 
physical feature of a given system, there are several argu- 
ments that advise caution in the interpretation of DFT 
calculations where this effect is found. The use of an 
approximate functional to model the exchange and cor- 
relation energy, such as LDA or GGA, generally produces 
severe and systematic errors in the values of (pp and </)„, 
which can be generally traced back to the well-known 
band-gap problemi^i^ This implies that finding a neg- 
ative value of either (pp or (p„ is unlikely to be a robust 
result of an LDA or GGA calculation. Furthermore, the 
total amount of spilled-out charge depends on the DFT 
values of (pp and (the more negative the SBH, the 
larger the number of states of the insulator that cross 
the Fermi level). This means that, in such a pathologi- 
cal regime, the error in (pp or (p„ will directly propagate 
to the charge density, and potentially affect a number of 
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fundamental ground-state properties of the interface. In 
order to avoid undesirable artifacts in the DFT results, 
it is therefore crucial to clearly identify whether this sce- 
nario applies to a given interface calculation. 

Such an analysis is not entirely straightforward, as 
physics governing the band alignment in a ferroelectric 
capacitor significantly departs from the well-established 
concepts of semiconductor physics. First, the imperfect 
screening at the electrode interface produces a poten- 
tial drop^^-^'^ that is roughly linear in the polarization 
P,^* and modifies the lineup between the bands of the 
insulator and the Fermi level of the metal<^ This phe- 
nomenon, central to the physics of ferroelectric capac- 
itors, has important implications for the stability of a 
monodomain polar statej^ and for devices based on the 
tunneling electroresistance efi^ect<^ Second, the residual 
"depolarizing" electric field produces a linearly increas- 
ing electrostatic potential in the film. This prevents a 
precise determination of the band lineup,— as a proper 
(and physically meaningful) definition of the latter re- 
quires a macroscopically constant reference energy in the 
insulating region. Third, the marked covalent character 
of bonding in perovskites produces non-trivial changes 
in the band structure of the insulator, depending on the 
magnitude of the polar distortion. This further compli- 
cates the extraction of an accurate band lineup by means 
of standard first-principles procedures, as the bulk refer- 
ence calculation needs to accurately match the electrical, 
in addition to the mechanical, boundary conditions of 
the film. Finally, and most importantly, one must keep 
in mind that all these new physical ingredients may co- 
exist with the more traditional features that are typical 
of metal/semiconductor interfaces, e.g. the phenomenon 
of metal-induced gap states (MIGS)<^ To guide future 
works in this field, and to build a firm theoretical basis 
for the interpretation of the experiments, it is becoming 
increasingly urgent to rationalize all these many compet- 
ing effects into a coherent picture, where the limitations 
of the current simulation methods can be clearly drawn. 

Here we develop a general and intuitive model of the 
band offset at a ferroelectric/metal interface, and its de- 
pendence on the polarization. We identify two quali- 
tatively distinct regimes, corresponding to (i) that of a 
normal Schottky alignment and (ii) that of a pathological 
Ohmic junction. We demonstrate the artifacts typically 
associated with (ii) by performing extensive calculations 
of technologically relevant ferroelectric/metal interfaces. 
We discuss the relevant literature works, pointing out 
those where our results suggest a revision of the currently 
accepted interpretation. We further identify a direct re- 
lationship between the pathological Ohmic regime and 
the physics of "electronic reconstruction"— at polar ox- 
ide interfaces such as LaAlOs/SrTiOa, and trace a viable 
route towards a unified description of these two phenom- 
ena. Finally, we discuss a number of viable methodolog- 
ical perspectives to overcome the limitations of DFT il- 
lustrated in this work. 

The paper is organized as follows: In Sec. [ll] 



we develop our theoretical model of the band off- 
set at a ferroelectric/metal interface, illustrating the 
main consequences of a "pathological" band align- 
ment. In Sec. mil we present a self-contained 
overview of the theoretical methods we use to de- 
tect such features in a first-principles calculation. In 
Sec. IIVI we present the results of our simulations for 
paraelectric capacitors, by comparing non-pathological 
(PbTiOa/SrRuOa and BaTiOa/SrRuOs) and pathologi- 
cal cases (KNbOa/SrRuOa and BaTiOa/Pt). In Sec. |V] 
we demonstrate that the two cases which we find to 
be non-pathological in the paraelectric configuration in- 
deed become pathological when the polarized ferroelec- 
tric state is fully relaxed. In Sec. lVII we discuss the impli- 
cations of this work with respect to the existing literature 
on the subject. Finally, in Sec. IVIII we present our con- 
clusions and outlook for future research. 



II. GENERAL THEORY OF THE BAND 
OFFSET 

A. Metal/semiconductor interfaces 

The Schottky barrier, a rectifying barrier for electrical 
conduction across a metal/semiconductor junction, is of 
vital importance for the operation of any modern elec- 
tronic device. For the case of an n-type semiconductor, 
the Schottky barrier height is the energy difference be- 
tween the conduction band minimum and the Fermi level 
across the interface, and we indicate it as The nature 
of the microscopic mechanisms governing the magnitude 
of (j)n has troubled scientists for several decades. In spite 
of the ongoing debates, it seems to be widely accepted 
now that, while bulk material properties certainly play a 
substantial role, 0„ is best understood as a genuine in- 
terface property. This is in agreement with the intuitive 
picture one gets from quantum mechanics: the charge 
rearrangement due to chemical bonding at the interface 
produces an interface dipole, and this will uniquely deter- 
mine the offset between the energy bands of the insulator 
and the Fermi level of the metal. 

To be more specific, it is useful to consider the elec- 
trostatic Hartree potential at the interface between two 
semi-infinite solids, 

^H(r) = / J^fr', (1) 

where p(r) is the total charge density (including elec- 
trons and nuclei). Vh is a rapidly varying function of the 
position, reflecting the underlying atomic structure. In 
order to filter out the large oscillations and preserve only 
those features that are relevant on a macroscopic scale, 
it is convenient to apply an averaging procedure i ^^i^^ 
This consists of (i) performing a global average of Vh (r) 
over planes parallel to the interface, and (ii) convoluting 
the resulting one-dimensional function with a Fourier fil- 
ter to suppress the high spatial frequency components. 
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FIG. 1: Schematic representation of the band offset at a 
metal/insulator junction, illustrating the main quantities dis- 
cussed in the text. 




FIG. 2: Schematic representation of a symmetric short- 
circuited ferroelectric capacitor in a polarized configuration 
within the imperfect-screening model, t is the thickness of 
the ferroelectric film. M and FE represent, respectively, the 
metal electrode and the ferroelectric film. Both materials are 
assumed to be separated by a vacuum layer of thickness Aoff. 
The thick solid line indicates the opposite of the electrostatic 
potential, ~V{z). 



(See Ref . ISOl for a detailed description of the method, 
and Ref. |3l| for an extensive review of its applications to 
SBH calculations.) After this "nanosmoothing"— proce- 
dure, the doubly-averaged Vniz) reduces to a step func- 
tion, from vifhich we can extract the electrostatic lineup 
termM^ 



diclcctricx 



ctalx 



(2) 



which includes all the physics of the interface dipole for- 
mation. [^Vjdioicctric^ ^y^otai^ ^j.g asymptotic val- 
ues of Vy{_{z) far from the interface.] To determine the 
band offsets from A(T^) it is then necessary to know how 
the bulk energy bands of the insulator and the Fermi 
level of the metal are related to their respective average 
electrostatic potential. In full generality, one can write 



j)p^EY~EY + ^{V), 
Ec-Ef + A{V). 



(3a) 
(3b) 



Ey, Ec and E-p are usually referred to as the band struc- 
ture termr^^ and are bulk properties of the two ma- 
terials. They are defined as the energy positions of the 
valence (Ey) and conduction (Eq) band edges of the in- 
sulator, and the Fermi level of the metal (Ep), all referred 
to the average (Vu) in the respective bulk (see Fig. [T]). 

In Sec. |TTT] we provide further details of the standard 
computational procedures used to calculate these quan- 
tities in practice. In the following Section we discuss how 
the above theory needs to be revised and extended in the 
case of metal/ferroelectric interfaces. 



B. Metal/ferroelectric interfaces 

Ferroelectric materials entail a new degree of freedom, 
the macroscopic polarization P, which is absent in the 
semiconductor case. It is natural then to expect that the 
above picture of the band offset at metal/insulator inter- 
faces may need to be extended to take this new variable 



into account. In the following, we discuss how P affects 
both the lineup and the band-structure terms in Eqs. ((5a|) 
and (l3bll. 



1. Lineup term 

We represent a simple ferroelectric material as a non- 
linear dielectric, which in bulk is characterized by an in- 
ternal energy Uh per unit cell of the form 



C/b(i?) = Ao + + AiD^ + 0{D^). 



(4) 



Here D is the electric displacement field, Aq is an arbi- 
trary reference energy, A2 is negative and the highest ex- 
pansion coefficient positive. (As we are concerned with 
the essentially one-dimensional case of a parallel-plate 
capacitor, we only consider the component of the D vec- 
tor that is normal to the interface plane, indicated as D 
henceforth.) The ^0,2,4,... coefficients implicitly contain 
all the complexity of the microscopic physics, and can 
be calculated from first principles using the methods of 
Ref. W. It follows from elementary electrostatics^^ that 
the internal electric field, £{D), is the derivative of U{D), 



1 dUi,{D) 
fl dD '■ 



(5) 



where Q. is the cell volume. 

The electrostatics of a parallel-plate capacitor config- 
uration can be well described^'2^'^ within the imperfect 
screening model, as sketched in Fig. [5] The iV-layer thick 
ferroelectric film can be thought of as separated from the 
ideal metal electrode by a thin layer of vacuum, of thick- 
ness Aoff. Aoff is an "effective screening length" that takes 
into account all the microscopic details of the interface 
dipole response to the polarization,— including electronic 
and chemical bonding effects;^ At the interface between 
the ferroelectric and the vacuum layer D must be pre- 
served. Therefore, an homogeneous electric field appears 
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Ds (C/m^) 


Aeff (A) 


A<^ (V) 


BaTiOs 


0.39 


0.20 


1.8 


PbTiOs 


0.75 


0.15 


2.6 



TABLE L Estimation of the change in the lineup term A(j) 
of typical ferroelectric capacitors upon polarization reversal. 
Ds is the bulk spontaneous polarization of the ferroelectric 
material. XcS were calculated in Ref. for capacitors with 
SrRuOa electrodes. 



inside the vacuum layer, of magnitude fvac — D/eg. Re- 
calling that the energy density of a static electric field £ 
in vacuum is u = — £>^/2eo, the energy of the 

iV-layer thick ferroelectric film can then be written as 

t/jvp) -iV(7bp) + 25Aeff^, (6) 

ZeQ 

where S is the surface cell area. [Note that two symmet- 
ric electrodes of equal Aog are considered in Eq. (|6]).] The 
second important consequence of a non-zero Aeff is that 
the lineup term, Eq. ([2]), now linearly depends on the ex- 
ternal parameter D, due to the additional potential drop 
at the interface, that can be computed as the product of 
the electric field within the vacuum layer times its width, 

A(\/)(i^) = A(V^)(0) + Aeff-. (7) 

Co 

[It is worth noting that, whenever £h{D) ^ 0, at the mi- 
croscopic level A{V){D) contains an intrinsic arbitrari- 
ness; furthermore, in such a case it is no longer justified 
to think in terms of an "isolated" interface between two 
semi-infinite solids. Techniques to deal with these issues 
in practical calculations are described in Ref. [20l .] 

To give a more quantitative flavor of the impact of this 
D-dependence in real systems, we can use the values of 
Aeff reported in the literature for PbTiOa/SrRuOs and 
BaTiOa/SrRuOs capacitors. Upon polarization reversal, 
the interface lineup term A{V) will undergo a variation 
corresponding to 

Acb^A{V)iDs)-A{V){-Ds)^2X,s — , (8) 

where Dg is the spontaneous polarization of the ferro- 
electric material (in the spontaneous configuration the 
internal electric field within the ferroelectric, 5b, vanishes 
and Dg equals the spontaneous polarization.) The values 
reported in Table U indicate that this effect can be rather 
large, of the order of 1-2 eV, even for ideal defect-free 
interfaces. 



2. Band- structure term 

The polar displacements in the ferroelectric film mod- 
ify not only the lineup term, but also the bulk band- 



structure term. This is most easily understood by re- 
calling the role played by covalent bonding in the ferro- 
electric instability of perovskite titanates. Hybridization 
effects between the cation 3d states and the oxygen 2p 
states are intimately linked to the off-centering of the Ti 
sublattice. This implies that the polar distortions can sig- 
nificantly modify both the conduction and valence band 
structure. For example, in both BaTiOa and PbTiOs 
the fundamentamental gap increases when going from the 
centrosymmetric cubic structure to the polar tetragonal 
phase. Using the arguments of Ref. [13, we can think of a 
continuous dependence of both Ey and Ec , respectively 
in Eq. pa)) and Eq. ([3b|. on the electric displacement D. 
The Fermi level Ep, of course, remains fixed as the elec- 
tric displacement does not affect the bulk of the metallic 
electrode. In summary, the general expression for the n- 
type SBH at a metal/ferroelectric interface (an analogous 
expression follows for the p-type one) is 

MD)^Ec{D)-Ef + A{V){D), (9) 

where at the lowest order Eq is quadratic in D (the lin- 
ear order is forbidden by symmetry), and in most cases of 
interest A{V){D) can be approximated by a linear func- 
tion as in Eq. ([7]). In the following, we shall elaborate on 
this expression and identify a new, qualitatively differ- 
ent regime, with important implications for the physics 
of the interface. 



C. Ferroelectric capacitors in a pathological regime 

Equation Q implies that (j3„{D) might become nega- 
tive for some values of D. From the point of view of first- 
principles calculations, already by looking at the values 
of TablcUwc can be reasonably sure that this will happen 
at the PbTiOa/SrRuOa interface: 2.6 eV is already larger 
than the LDA gap of PbTiOa in the ferroelectric phase 
(^^2.0 eV). This possibility has been almost systemati- 
cally overlooked in the literature. As this is a central 
point of this work, we shall illustrate in detail the conse- 
quences of such a regime, and explain why we regard it as 
"pathological" . We discuss in the following two possible 
occurrences of this scenario: (i) (/)„ is negative already 
in the paraelectric configuration at D = and (ii) is 
positive at -D = but becomes negative at some value of 
\D\<Ds. 



1. The centrosymmetric case 

We start with a capacitor in the reference paraelec- 
tric structure with two symmetric electrodes, and we hy- 
pothesize that, for whatever physical reason, the interface 
dipole that forms between the metal and the film leads to 
a negative 0„. (Similar arguments apply to the case, not 
explicitly discussed here, of a negative 0p.) As the quan- 
tum states of the conduction band of the film lie at lower 
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energy than the Fermi level of the metal, the former will 
be filled up to Ep, leading to a nonzero free charge den- 
sity, pfrcG, in the film. Neglecting quantum confinement 
effects, we can use the Thomas-Fermi model and treat 
the free charge distribution as macroscopically uniform. 
Within this approximation, ptroc is exactly given in terms 
of 4>n and the electronic density of states of the bulk insu- 
lator, Ph{E), in a vicinity of the conduction band bottom. 



(a) 



Pfr, 



Ec-e4 



Ec 



Ph{E)dE. 



(10) 



This additional charge density, superimposed on an oth- 
erwise charge-neutral insulating film, will produce a 
strong electrostatic perturbation in the system. For ex- 
ample, if such a charge rearrangement occurred in vac- 
uum, the Poisson equation 



(PV{z) 



Pfrcc 



would imply a parabolic potential of the form 

T// \ Pfrec 2 

y{z) = -^r—z . 
Zeo 



(11) 



(12) 



(We assume that z = corresponds to the center of the 
ferroelectric film.) Throughout this work, we shall as- 
sume that the interface is oriented along the z axis, and 
each material is periodic in the plane parallel to the in- 
terface, referred to as the {x, y) plane. As typical ferro- 
electric materials are exceptionally good dielectrics, in a 
first approximation we can assume that V{z) will be per- 
fectly screened by the polar displacements of the lattice. 
However, this does not mean that electrostatics has no 
consequences - quite the contrary. Macroscopic Maxwell 
equations in materials indeed dictate that 



dD{z) 
dz 



(13) 



Hence, if we assume perfect bulk screening, we have 
£{z) = 0, D{z) = eo£{z) + P{z) = P{z) and, after inte- 
grating Eq. (IT3l) . P{z) = pfieeZ- So, since the sign of the 
electronic charge and phee are negative within our con- 
vention, we have a non-uniform and linearly decreasing 
polarization in the ferroelectric film [see Fig.[31[d)]. This 
means that, at the film boundaries (z = ±i/2, where t 
is the thickness), the local electric displacement has now 
opposite values, proportional to the total amount of free 
charge that was transferred. 



D( 



:Pfi< 



(14) 



Of course, the band offset at the interface depends on the 
local value of D in the film region adjacent to the inter- 
face, so 4>n will be consequently shifted in energy accord- 
ing to Eq. (jni). We can expect that for small D values 
the (quadratic) polarization effects on the band structure 



(c) 



(b) 

-I- 
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1/2 z 




FIG. 3: Schematic representation of the effect of free-charge 
redistribution onto the band diagram of a paraelectric capac- 
itor with a negative (a) band alignment under perfect 
interface screening (i. e. when pfree vanishes), and (b) after 
charge spill-out and electrostatic re-equilibration. The corre- 
sponding profiles of the electric displacement field within the 
ferroelectric films are displayed in panels (c) and (d), obtained 
after integrating Eq. p3|) . 



will be less important than the (linear) dependence of the 
lineup term on D. (Note that the presence of additional 
charge in the conduction band might also alter the band- 
structure term, e.g. through on-site Coulomb repulsions 
or other exchange and correlation effects; in the limit of 
weak correlations we expect these to be even smaller and 
essentially irrelevant for this discussion.) Therefore, we 
approximate Eq. ^ with Eq. ([7]), and write 



Co 



2en 



(15) 



[The minus sign comes from the fact that at the z < 
interface, which is the one for which Eq. ^ is valid within 
our conventions, D is negative.] In turn, the new 0„ will 
modify pfroc through Eq. (jlOl) . For some value of </>„, 
Eq. ([T0| and Eq. (|T5|) will be mutually self-consistent 
and the system will reach electrostatic equilibrium. This 
can be expressed through an integral equation where we 
have eliminated phcc 



Ph[E)dE = ■ 



Ec 



(16) 



To qualitatively appreciate the physical implications 
of this expression, we can explicitly solve it by using a 
constant ph{E) = a. (Note that this assumption is not 
completely unrealistic as the t2g bands forming the bot- 
tom of the conduction band in many ferroelectric per- 
ovskites have a marked 2D character; in other words, the 
in-plane effective mass mjj is much smaller than the out- 
of-plane one, m*^. Within the approximation m*^ = oo, 
the constant density of states of the six-fold degenerate, 
free-electron- like 2D band is uniquely determined by to|.) 
This leads to 



(17) 
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and with a few rearrangements to 

C-rAcffa + 1 

where C = e^/2eo is a constant, and a = a/Q. is the 
density of states per unit energy and volume of the bulk. 
In spite of the drastic simplifications, Eq. (|18p already 
contains most of the relevant ingredients for our analysis. 
A few notable ones are missing - we shall come back to 
those in Sections llI C 2l and lII C 31 Before going into more 
detailed considerations, however, it is important to spell 
out the direct implications of Eq. ([T8| , which we shall be 
concerned with in the following. 

First, note that all quantities appearing at the denom- 
inator at the right-hand side of Eq. ([T8| are positive. 
This means that will be negative, and will satisfy 
< 0n < 0. The lower limit corresponds to the per- 
fect interface screening case, AcS = 0. The upper limit 
corresponds to no screening, Acs — > oo. The situation 
is schematically represented in Fig. [SJa) and Fig. [Hb). 
Given a negative 0° [Fig.EIa)], the charge redistribution 
will induce an upward energy shift of the conduction band 
minimum (CBM), bringing (/)„ closer to the Fermi level 
[Fig. IHl^b)]. Second, in the limit of i — >■ cx) (infinite thick- 
ness) (pn will tend to zero from below as oc —1/t. This 
means that the self-consistent band offset 0„ is not de- 
termined by the local physical properties of the junction, 
i.e. it is no longer an interface property - the spilled-out 
charge will redistribute over the whole film thickness as 
t is varied. Third, the density of states of the conduction 
band, represented in Eq. ((T5| by the parameter a, will 
also affect the value of 0„: the larger a, the strongest 
the reduction in (/)„ upon charge spill-out and electro- 
static re-equilibration. (To avoid confusion, note that in 
the above paragraphs, we used the word "screening" in 
two different contexts. By "perfect bulk screening" we 
mean £h{D) = 0. By "perfect interface screening" we 
mean Aoff =0.) 

We can attempt a semiquantitative assessment of 
Eq. in a representative capacitor of thickness t = 50 
A (comparable to those that are typically simulated 
within DFT). In atomic units, we use AeS — 0.3 (of 
the order of the values reported in Table U, C = 27r, 
and a — 0.05 (appropriate for the conduction band of 
SrTiOs, a prototypical perovskite material, with a calcu- 
lated mji = 0.77 and f2 = 385 a.u.). We obtain 

K - (19) 

This implies that the effect is quite strong - even if 0" 
is a rather large negative value (e.g. of the order of -1 
eV) , in most practical cases the conduction charge redis- 
tribution will reduce it to a value that lies just below the 
Fermi level. Most importantly, this implies that, when 
(/)° < 0, the physical parameters, (/)° and Aeff, governing 
the band offset at the interface are neither accessible in 
a simulation, nor are they directly measurable in an ex- 
periment - only (/)„ might be. Note, however, that the 



"self-consistent" (/>„ value is generally not a well-defined 
physical quantity ~ this is only true within the many 
approximations used in the above derivations. In partic- 
ular, we have neglected band-bending effects: in general 
the electrostatic potential will be non-uniform in the film 
(see Sec. Ill C 31) and 0„ will be a function of the distance 
from the interface. But even if we put this caveat aside 
for a moment, the reader should keep in mind that (f>n 
is determined here by space- charge effects through sev- 
eral independent contributions. Furthermore, the film 
is no longer insulating but becomes a metal. This is a 
substantial, qualitative departure from the physical con- 
cepts that were developed in the context of semiconduc- 
tor/metal interfaces, and that led to the consensus un- 
derstanding of (f>n as a genuine interface property. 

Given this situation, one needs to revisit the very foun- 
dations of the methodological ab-initio approaches that 
have been used with great success in the past to com- 
pute Schottky barrier heights. This success has critically 
relied on a key observation: the interface dipole, that 
one identifies with the lineup term Eq. is a ground- 
state property, i.e. is not directly affected by the well- 
known limitations of the Kohn-Sham eigenvalue spec- 
trum. This is excellent news: one can efficiently (and 
accurately) calculate A{V) within DFT, and combine it 
with a band-structure term (Ey or Eq) calculated at a 
higher level of theory (e.g. GW); within this formally 
sound procedure, theoretical calculations have shown re- 
markable agreement with the experimental observations 
in the past. 

In the spill-out regime (i.e. 0° < 0) described in this 
Section the above key observation no longer holds - the 
erroneous DFT value of 0^ plays a direct and dominant 
role in the interface dipole formation, as is apparent from 
Eq. Furthermore, as 0^ is systematically underesti- 
mated within LDA or GGA, there is the concrete possi- 
bility that the spill-out regime itself (0° < 0) might be an 
artifact of the band-gap problem. Thus, the ground-state 
properties of the system found in a simulation might be 
qualitatively wrong due to this issue, in loose analogy to, 
e.g., the erroneous LDA prediction of metallicity in many 
transition metal compounds. It goes without saying that 
the results of a simulation where significant spill-out of 
charge is found because of the mechanism described in 
this Section should be regarded with great suspicion. 



2. The broken-symmetry case 

Even if the band aligmnent is Schottky-like in the ref- 
erence paraelectric structure of the capacitor, Eq. ^ en- 
tails the possibility that it might become pathological in 
the ferroelectric regime (i.e. when the polar instability is 
allowed to fully relax). Unfortunately, for this case many 
of the simplifying assumptions used above are no longer 
valid, and for a detailed description one would need to 
take into account the more refined physical ingredients 
discussed in Sec. Ill C 31 At the qualitative level, how- 
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FIG. 4: (Color online) (a) Paraelectric capacitor with a 
Schottky-like band alignment in the paraelectric structure, 
(b) When the polar instability sets in, the band alignment 
becomes pathological, the conduction band is locally popu- 
lated (red shaded area) and the film becomes partially metal- 
lic (light shaded area bounded by the dashed line). 



ever, we can already draw some important conclusions, 
as we shall briefly illustrate in the following. 

Eq. ([9|) predicts that, if 0^ is positive and the capac- 
itor is compositionally symmetric [as in Fig. UJJa)], at fi- 
nite D at most one of the two opposite interfaces will 
have a negative 0„. This implies that only part of the 
ferroelectric film, i.e. the region adjacent to this "patho- 
logical" interface, will become metallic, while the rest of 
the film will stay insulating [Fig. IH^b)]. (To understand 
this point, note that in contrast with the previous case 
one has now a finite "depolarizing" electric field in the 
insulating part of the capacitor. This wedge-like poten- 
tial will keep the conduction electrons electrostatically 
confined to the pathological side.) In the insulating re- 
gion, the polarization will be macroscopically constant, 
as in a well-behaved capacitor [recall Eq. ([TB])]. Accord- 
ing to the same Eq. (|13|) . D{z) [and hence P{z)] will be 
non-homogeneous, with a negative slope, in the metallic 
region. 

In this context it is worth pointing out an important 
physical consequence of such a peculiar electronic ground 
state. This concerns the response of the capacitor to an 
applied bias potential. In well-behaved cases, the po- 
larization of the capacitor will respond uniformly to a 
bias, i.e. all the perovskite cells up to the electrode in- 
terface will undergo roughly the same polar distortion. 
In the present "ferroelectric-pathological" regime, part 
of the ferroelectric film has become metallic, i.e. the 
metal /insulator interface has moved to a place that lies 
somewhere in the film. This means that, if one tries to 
switch the device with a potential, the electric field won't 
affect the dipoles that lie closest to the pathological in- 
terface - they will be screened by the spilled-out free 
charge. A consequence is that the dipoles near a patho- 
logical interface will appear as if they were pinned to a 
fixed distortion, that is almost insensitive to the elec- 
trical boundary conditions. This pinning phenomenon 
has been studied in earlier theoretical works, and was as- 
cribed to chemical bonding effects. In Sec. |V] we shall 
substantiate with practical examples that "dipole pin- 
ning" is instead a direct consequence of the problematic 
band-alignment regime described here. In Sec. I VII we 
shall come back to this point and put it in the context of 
the relevant literature. 



3. Towards a quantum model 

In order to draw a closer connection between the 
semiclassical arguments of the previous sections and 
the quantum-mechanical results that we present in Sec- 
tions IIVI and |Vl we briefiy discuss here how to improve 
our physical understanding of the charge spill-out process 
by lifting some of the simplifying approximations used so 
far. As a detailed treatment goes beyond the scope of 
the present work, we shall limit ourselves to qualitative 
considerations. 

The most drastic approximation of our model appears 
to be the assumption of perfect dielectric screening within 
the ferroelectric material, where the spill-out charge is 
perfectly compensated by the polar displacements of the 
lattice. This implies that the electric field in the film 
vanishes, and the excess conduction charge can spread 
itself spatially at essentially no cost. In this scenario, the 
macroscopically uniform distribution of /Oftoe postulated 
in Sec. Ill C II appears very reasonable. In reality, the 
internal £ field in the bulk ferroelectric material does 
not vanish, but is a non- linear function of _D, that can be 
written by combining Eq. ^ and Eq. ([S]), 

£h{D) ^ ^ {2A2D + AAiD^) (20) 

Of course, solving for the self-consistent pheeiz) in a 
non-linear medium would require a numerical treatment. 
Still, we can gain some insight about qualitative trends by 
starting, for example, from the linearly decreasing D{z) 
found in the D = case of Sec. Ill C II Using Eq. ^ we 
can write £{z) — £b[D{z)]. The electrostatic potential is 
then given by integrating £{z). This essentially leads to 

V]i{z) = i7b[-D(z)]/(5oj where U\, is the internal electro- 
static energy of Eq. Q , and Qo is a (positive) constant 
with the dimension of a charge. This means that the 

spatial variation in V-a_{z) reflects the energy landscape 

of the bulk material: V-niz) will be a double- well poten- 
tial in a ferroelectric material {A2 < 0), and a single- well 
potential in a paraelectric material {A2 > 0). Remark- 
ably, the double-well potential accounts for the possi- 
bility of free-charge accumulation in the middle of the 
centrosymmetric film [Fig.[5ja)], which would produce a 
head-to- head domain wall in the polarization P{z). Con- 
versely, for a paraelectric material one would expect the 
free charge to be (more or less loosely) bound to the in- 
terface, and have a minimum in the middle of the film 
[Fig. Eljb)]. Of course, these considerations are valid for 
a centrosymmetric capacitor, and are presented just to 
give the reader an idea of the physics - in the ferroelec- 
tric case, more complex patterns can occur and exploring 
them all would require an in-depth study that is beyond 
the scope of this manuscript. 

A second important approximation is the neglect of (i) 
quantum confinement effects beyond the simple Thomas- 
Fermi filling of the bulk-like density of states and (ii) 
the band-structure changes due to the polar distortions, 
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FIG. 5: Schematic representation of the efTects of dielec- 
tric nonlinearity on the band diagram of a centrosymmetric 
capacitor. The effective potential felt by the conduction elec- 
trons is ~Vh{z) (see text), (a): Ferroelectric material, (b): 
Paraelectric material. 



which we briefly mentioned in Sec. Ill B 21 These will 
further modify the equilibrium distribution of the free 
charge, and we expect them to be important to gain a 
truly microscopic understanding of the system, although 
not essential for the points of this work. Remarkably, a 
promising model taking all these ingredients into account 
(dielectric non-linearity and band-structure effects) was 
recently proposed in the context of the (at first sight un- 
related) LaAlOs/SrTiOa interface.'^* This indicates that 
the physics of a ferroelectric capacitor in the patholog- 
ical band-alignment regime described here is essentially 
analogous to that of the "electronic reconstruction"— in 
oxide superlattices. Further work to explore these inter- 
esting analogies is under way. 



tem. Note that this feature has been ascribed in earlier 
works to phenomena of completely different physical ori- 
gin. We shall devote special attention in Sections ITVl and 
IVl to demonstrating the intimate relationship between 
Pfrco and spatial variations in P. 



III. METHODS 

In this Section we spell out the practical techniques 
that we use to extract the SBH from first-principles cal- 
culations, the operational definitions of free charge and 
bound charge, and the methods we use to control the 
electrical boundary conditions in supercell calculations. 
We also summarize the other relevant computational pa- 
rameters used in Sections IIVI and IVl 



A. Schottky barrier estimations 

First, we briefly review the methods that were 
used in earlier works to compute Schottky barriers at 
metal/semiconductor interfaces, pointing out advantages 
and limitations of each of them. Then, we illustrate po- 
tential complications that might arise, with special focus 
on ferroelectric oxide systems and the issues discussed in 
Sec. EH 



D. Implications for the analysis of the ab-initio 
results 

The above derivations show that there are two 
qualitatively dissimilar regimes in the physics of a 
metal/insulator interface, Ohmic-like and Schottky-like. 
During the derivation, we have evidenced some distinct 
physical features that we expect to be intimately associ- 
ated with the "pathological" Ohmic case. As these are of 
central importance to help distinguish one scenario from 
the other, we shall briefly summarize them in the follow- 
ing, mentioning also how each of these "alarm flags" can 
be detected in a first-principles simulation. 

First, even after the electron re-equilibration takes 
place, the band edges cross the Fermi level of the metal, 
i.e. the apparent Schottky barrier is negative. Therefore, 
the analysis of the local electronic structure and of the 
SBH appears to be the primary tool to identify a patho- 
logical case. However, as the "self-consistent" </)„ tends 
to stay very close to the Fermi level, this analysis should 
be performed with unusual accuracy - techniques to do 
this will be discussed in Sec. IIII Al 

Second, the presence of a substantial density of free 
charge populating the conduction band of the insula- 
tor is another important consequence of the pathologi- 
cal regime. In Sec. IIII BT] we illustrate how to rigorously 
define pfroo in a ferroelectric heterostructure. 

Finally, a remarkable consequence of charge spill out is 
the presence of an inhomogeneous polarization in the sys- 



1. From the local density of states 

In order to calculate the band offset at a 
metal/insulator interface, one needs to identify the lo- 
cation of the band edges deep in the insulating region, 
with the Fermi level of the metal taken as a reference. 
To that end, it has become common practice^S, to define 
a spatially-resolved density of states, 

p(»,S)=V/ dk|(i|V;„k)l''^(^-£^nk), (21) 
n -'BZ 

where \i) is a normalized function, localized in space 
around the region of interest. When \i) = |r) is an 
eigenstate of the position operator, the resulting p(r, E) 
is commonly known as local density of states (LDOS). 
Conversely, when \i) = \4>nim) is an atomic orbital of 
specified quantum numbers {nj,m), we call it instead 
projected density of states (PDOS). The integral is per- 
formed over the first Brillouin zone (BZ) of the supercell 
and the sum runs over all the bands n. Enk stands for 
the eigenvalue of the electronic wave function V'nk- 

The LDOS defined in Eq. that depends on the 

position in real space as well as on the energy, gives a 
very intuitive picture of the band offset: "sufficiently far" 
away from the interface, the LDOS converges to a bulk- 
like curved and in principle the location of the band 
edges (and hence the SBH) can be directly extracted 
by visual inspection. However, several approximations 
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are used in practice to make the calculation tractable, 
and these can introduce significant deviations in the SBH 
computed by means of either the LDOS or PDOS. First, 
all studies are done on a finite supercell, usually with 
a symmetric capacitor geometry. This implies that the 
LDOS of the most dispersive bands will be altered by 
quantum confinement effects, which might produce a spu- 
rious gap opening. Also, the LDOS associated to the 
evanescent metal-induced gap states (MIGS) might be 
still important at the center of an insulating film that is 
not thick enough, thus preventing an accurate identifica- 
tion of the band edge. Second, a discrete fc-point mesh is 
used instead of the continuous one implicitly assumed in 
Eq. ([2T|) . Such a fc-point mesh is generally optimized for 
efficiency, which means that high-symmetry (HS) points 
are often excluded. '^^ As the edges of the valence and 
conduction band manifolds are usually located at the HS 
points, "^^ estimating those features from the calculated 
LDOS might lead to substantial inaccuracies. For ma- 
terials that display a very dispersive band structure (see 
e.g. Ref. HI) it is not unusual to have deviations of the 
order of several tenths of an eV. Third, a fictitious elec- 
tronic temperature (or Fermi surface smearing) is com- 
monly used, in order to alleviate the errors introduced 
by the fc-mesh discretization. This implies that the Dirac 
delta function in Eq. (|2ip needs to be replaced by a nor- 
malized smearing function (e.g. a Gaussian) with finite 
width. This is a again potential source of inaccuracies, 
because the apparent edges of the smeared LDOS/PDOS 
actually might not correspond to the physical band edges 
but to the (artificial) tail of the smearing function used. 

Summarizing the above, we get to the following oper- 
ational definition of the smeared LDOS, 

p(r,S) = ^Wk|V^„k(r)|'5(^-^nk), (22) 

nk 

where the BZ integral has been replaced with a sum 
over a discrete set of special points k with corresponding 
weights ifkj and the Dirac delta has been replaced with 
a smearing function g. As will become clear shortly (a 
detailed analysis is provided in Appendix [B|) . it is very 
important to use in Eq. (|22p a g function that is minus 
the analytical derivative of the occupation function used 
in the actual calculations. The Gaussian smearing (G) 
and the Fermi-Dirac (FD) smearing are by far the most 
popular choices. These correspond to the following defi- 
nitions of g, 

g^{x)^^e~-'/-\ (23a) 

where a is the smearing energy used during self- 
consistent minimization of the electronic ground state. 



2. From the electrostatic potential 

To work around these difficulties, it is in most cases 
preferrable to avoid the direct estimation of the SBH 
based on the LDOS/PDOS, and use instead the indirect 
procedure, based on the nanosmoothed electrostatic po- 
tential, Fh, described in Sec. IH Al The interface lineup 
term, A(V^), generally (a notable exception is the patho- 
logical spill-out regime described in Sec. [H]- for further 
details see Sec. IHI A4|) converges much faster than the 
LDOS/PDOS with respect to all the computational pa- 
rameters described above (slab thickness, fc-mesh, Fermi 
surface smearing). The band-structure terms, Ey and 
i?c, can be then accurately and economically evaluated 
in the bulk, without the complications inherent to MIGS 
and quantum confinement effects. While this is in princi- 
ple a very convenient and robust methodological frame- 
work it is, however, also prone to systematic errors. In 
particular, great care must be used when performing the 
reference bulk calculations. In the vast majority of cases 
these must not be performed on the equilibrium structure 
of the bulk solid, but will be constructed to accurately 
match (i) the mechanical and (ii) the electrical boundary 
conditions of the insulating film in the supercell. The 
issue (i) is well known: in a coherent heterostructure the 
insulating film is strained to match the substrate lat- 
tice parameter, and for consistency the "bulk" calcula- 
tion should be performed at the same in-plane strain. 
(The dependence of the band-structure term on the lat- 
tice strain is well known in the literature, and referred 
to as "deformation potentials''.'^^) Issue (ii) concerns fer- 
roelectric systems, and is therefore not widely appreci- 
ated within the semiconductor community. Whenever 
the symmetry of the capacitor is broken and there is a 
net macroscopic polarization in the ferroelectric film, the 
structural distortions may alter the band structure signif- 
icantly, often more than purely elastic effects do.— Note 
that in most capacitor calculations the film is only par- 
tially polarized (i.e. it has neither the centrosymmetric 
non-polar structure, nor the fully polarized ferroelectric 
structure because of the depolarizing effects described 
in Sec. Ill Bp . The "bulk" reference calculation should 
then accurately match the polar distortions of the film^ 
extracted in a region where the interface-related short- 
range perturbations have healed into a regular pattern. 

3. The "best of both worlds" 

In order to minimize the drawbacks associated with ei- 
ther of the two methods described above, we find it very 
convenient to combine them in the following procedure. 
First, we compute the LDOS in the supercell at an atomic 
site (or layer) located far away from the interfaces, where 
the relaxed atomic structure has converged into a regular 
pattern. Second, we extract the relaxed atomic coordi- 
nates from the same region of the supercell, and build 
a periodic bulk calculation based on them, by preserv- 
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ing identical structural distortions and strains, and by 
using an equivalent fc-mesh. (An approximation is made 
here, since the periodic bulk simulation is carried out at 
zero macroscopic field while the LDOS in the supercell 
might contain the effects of a non-zero depolarizing field. 
The problem of computing the bulk layer-by-layer LDOS 
under a finite electric field remains an open question.) 
Third, we extract the LDOS from the bulk at the same 
atomic site or layer; we construct the bulk LDOS us- 
ing Eq. (j22p and an identical g function to that used in 
the supercell. Finally, we superimpose the bulk LDOS 
on the supercell LDOS at each layer j; we align them 
by matching the sharp peaks of a selected deep semi- 
core band, which are located at energies i?|c''°"^°''(j) ^'^d 
The deep semicore states are insensitive to the 
chemical environment and have negligible band disper- 
sion; this means that they provide an excellent, spatially 
localized reference energy for the estimation of the lineup 
term. 

At this point, we look at either LDOS curve in the 
vicinity of the Fermi level. If it is non-zero we are proba- 
bly facing a pathological spill-out case (see the following 
Section). If it is zero, then we can go one step further 
and accurately estimate the positions of the local band 
edges. To this end, we compute from the bulk calcula- 
tion £;g""^ and £;^""^, together with E^^^^. (A further 
non-selfconsistent run might be needed if the original fc- 
mesh did not include the HS fc-points where the band 
edges are located.) Finally, assuming that £'|^P°'''^°''(j) 
are all referred to an energy zero corresponding to the 
self-consistent Fermi level of the supercell, we define the 
local position of the band edges as 

i?cT"™^j) = - £^'c"" + (24) 

This procedure avoids the (often inaccurate) estimate of 
the band edges based on the tails of the smeared LDOS, 
and at the same time preserves the advantages of the 
"lineup -I- band structure" technique. In principle, the 
latter method should accurately match the results of 
Eq. (|24p . except for quantum confinement effects in the 
metallic slab used to represent the semi-infinite electrode, 
as discussed in Ref. |40. 

Note that this technique is not only useful to detect 
pathological band alignments and extract accurate band 
offsets in the non-pathological cases. Given that we are 
superimposing two LDOS calculated with identical com- 
putational parameters and structures, their direct com- 
parison can be very insightful. Most importantly, one 
expects all the features to closely match unless there are 
MIGS or confinement effects. Therefore, one has also a 
powerful tool to directly assess the impact of the latter 
physical ingredients in the supercell electronic structure. 
This procedure, therefore, yields far more physical infor- 
mation than the separate use of either the PDOS/LDOS 
or the nanosmoothing method. 



4- The pathological regime 

In the pathological regime described in Sec. [Ill many 
of the conditions that formally justify application of the 
above methods to the estimation of the SBH break down. 
First, the presence of a non- uniform electric displace- 
ment D{z) implies that the polar distortions are also non- 
uniform, and they may not converge to a regular bulk-like 
pattern anywhere in the film. Second, electrostatic and 
exchange and correlation effects due to the partial filling 
of the conduction band imply that the band structure 
may significantly depart from what one computes in the 
insulating bulk (note that this is distinct from the effect 
of the structural distortions discussed in the previous Sec- 
tion) . Third, the usual assumption of fast convergence of 
the interface dipole with respect to slab thickness, fc-mesh 
resolution and smearing energy also breaks down, as the 
conduction band DOS (which converges slowly with re- 
spect to these parameters) is now directly involved in the 
electrostatic re-equilibration process. Based on this, the 
reader should keep in mind that there is an intrinsic arbi- 
trariness, of physical more than methodological nature, 
in the definition of the band edges in spill-out cases. This 
arbitrariness reflects itself in the fact, already pointed out 
in Sec. [TTl that the band alignment at a pathological in- 
terface is no longer a well-defined interface property, nor 
is it directly measurable in an experiment. The position 
of the bands is essentially the result of a complex elec- 
tron redistribution process that may occur on a scale that 
is almost macroscopic, and is driven by different factors 
than those usually involved in the SBH formation. 

Of course, by using all the precautions that are valid at 
well-behaved interfaces one might still gain some qualita- 
tive insight into the local electronic properties of the sys- 
tem. However, the data must be interpreted with some 
caution, and it is most appropriate to combine the anal- 
ysis with other post-processing tools before drawing any 
conclusion. We shall discuss some of these further anal- 
ysis tools in the following Sections. 



B. Electrical analysis of the charge spill-out 

In this Section we introduce the methodological tools 
that we use to analyze in practice the spill-out regime, 
in light of the theory developed in Sec. HT] In partic- 
ular, we illustrate how to rigorously define the "local 
electric displacement" D{z) and the "conduction charge" 
Pfioo- To evaluate the former, we discuss two approaches. 
The first one is based on a Wannier decomposition of 
the bound charges. The second one is an approximate 
formula in terms of the ionic distortions and the Born 
effective charges (BEG). This simplified formula is very 
practical for a quick analysis, but is generally affected 
by systematic errors. We address this issue by proposing 
a simple correction that significantly improves the accu- 
racy of the BEG estimate. 
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1. Definition of bound charge and conduction charge 

In a typical metal, it is difficult to rigorously identify 
conduction electrons and bound charges, as usually the 
respective energy bands intersect each other in at least 
some regions of the Brillouin zone. (This is true, for ex- 
ample, in all transition metals, where the delocalizcd sp 
bands cross the more localized d bands.) By contrast, 
in all perovskite materials considered here, even upon 
charge spill out and metallization, a well-defined energy 
gap persists between the bound electrons and the par- 
tially filled conduction bands. Therefore, it is straight- 
forward to separate the two types of charge densities, 
free and bound, simply by integrating the local density 
of states, defined in Eq. ([22]). over two distinct energy 
windows. For example, for the conduction charge phco 
we have 

Pfrco(r) = / p{T,E)dE= ^ Wk/„k |'0nk(r)|^ , 

E„k>Eo 

(25) 

where Eq is an energy corresponding to the center of 
the gap between valence and conduction band, p is the 
smeared density of states of Eq. ([^ , /„k are the occupa- 
tion numbers and the sum is restricted to the states with 
eigenvalue i?„k higher than Eq. [Note that Eq. (^5)) only 
holds if the (7-smearing of p is compatible with the def- 
inition of /„k, see the last paragraph of Sec. IIII A II and 
Appendix [BJ] Since we are working with layered systems 
that are perfectly periodic in plane, we will be mostly 
concerned with the planar average of pfree, 

Pfroo(2) = ^ phcc{r)dxdy, (26) 

where S is the area of the interface unit cell. In some 
cases, it is also useful to consider the nanosmoothed 
functiouj^S which we indicate by a double bar symbol, 

Pfree(^)- 

Concerning the bound charges, we shall approximate 
the local electric displacement D{z) with the local polar- 
ization P{z). This is an excellent approximation in many 
ferroelectric materials, where P is of the order of 0.1-1 
C/m~^ and D — P ^ eqE is typically much smaller than 
10"'^ C/m^^. (The largest electric fields £ that can be 
applied without dielectric breakdown^ ^ are of the order 
of 0.1 GV/m.) Thus, assuming D{z) ~ P{z) entails er- 
rors of 1% or less, which we consider negligible for the 
purposes of our discussion. Techniques to extract P{z) 
from a supercell calculation are described in the following 
sections. 



2. Local polarization via Wannier functions 

A very useful tool to describe the local polariza- 
tion properties of layered oxide superlattice are the 



"layer polarizations" introduced by Wu et al^ First, 
we transform the electronic ground state into a set of 
"hermaphrodite" Wannier orbitals^i^ by means of the 
parallel-transport''^ procedure. Note that we restrict the 
parallel-transport procedure only to the orbitals that we 
consider as "bound charge" , i.e. those with an energy 
eigenvalue lower than Eq. Then, we group the Wannier 
centers and the ion cores into individual oxide layers, and 
define the dipole density of layer j as 

where Za is now the bare valence charge of the atom a, 
whose position along z is Raz, and Zi is the location of 
the Wannier orbital i. 

Note that individual oxide layers in II-IV perovskites 
like BaTiOa or PbTiOa arc charge-neutral and the pj 
are well-defined; however, in I-V perovskites like KNbOs, 
individual layers are charged, and the pj become mean- 
ingless as they are origin dependent. To circumvent this 
problem, one can either combine the layers two by two as 
was done in Ref. |45| . or perform some averaging with the 
neighboring layers, as for example in Ref. 43- It is im- 
portant to keep in mind that, depending on the specific 
averaging procedure, one might end up with the formal 
or with the effective local polarization;^® in this work we 
find it more convenient to work with the latter. As we 
do not need, for the purpose of our discussion, to resolve 
P into contributions from individual AO and BO2 ox- 
ide layers, at variance with Ref. i43| we perform a simple 
average 

1 1 1 

We then define the local polarization by scaling this sur- 
face dipole density by the average out-of-plane lattice pa- 
rameter, c, of the oxide film, and by taking into account 
that every individual oxide layer occupies only half the 
cell. We thus define the local polarization as 

Pj = ^Pj- (29) 

The local polarization Pj is, of course, a discrete set of 
values, but we can think of it as a continuous function 
of the z coordinate, P{z), which is sampled at the oxide 
plane locations. In the remainder of this work, we will 
write Pj or P{z) depending on the context, but the reader 
should bear in mind that these two notations refer to the 
same object. 

3. Approximate formula via Born effective charges 

While the above definition of Pj in terms of Wannier 
functions is accurate and rigorous, it is not immediately 
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available in most electronic structure codes. An approxi- 
mate estimate of the local polarization can be simply in- 
ferred from the bulk Born effective charges and the lo- 
cal atomic displacements. Analogously to the above for- 
mulation, we can write the Z*-based approximate layer 
dipole density, pj, as 

pf = ^Y.K^o.z, (30) 

where Z* is now the bulk Born effective charge associ- 
ated with the atom a. Again, are ill-defined in per- 
ovskite materials, as typically individual oxide layers do 
not satisfy the acoustic sum rule separately. To address 
this issue, we perform an analogous averaging procedure 
and define 

Pj = -^Pj-i + ^Pj + -^Pj+i- (31) 

The approximate local polarization then immediately 
follows, 

= -pI (32) 

Such an approximation provides an exact estimate, in 
the linear limit, of the polarization induced by a small 
polar distortion under short-circuit electrical boundary 
conditions, i.e. assuming that the macroscopic electric 
field vanishes throughout the structural transformation. 
Neither of these conditions is respected in a ferroelec- 
tric capacitor, where the polar distortion is generally 
large (close to the spontaneous polarization of the fer- 
roelectric insulator), and where there is generally an im- 
perfect screening regime, with a macroscopic "depolariz- 
ing field" .— We investigate both issues in the Appendix, 
where we find that a simple scaling factor corrects, to 
a large extent, the discrepancy between Pj and . In 

particular, we write the "corrected" P^ as 

Pf =(l + pf^ (33) 

V XION/ 

where X(X) and Xion are, respectively, the electronic and 
ionic susceptibilities of the bulk material in the cen- 
trosymmetric reference structure, calculated at the same 
in-plane strain as the capacitor heterostructure. Note 
that for a ferroelectric material in the centrosymmetric 
reference structure, xion is negative, which is a conse- 
quence of the polar unstable mode in the phonon spec- 
trum. This means that the scaling factor will be smaller 
than 1 (^0.9 for the materials considered in this work). 
Practical methods to calculate X(xi and Xion are reported 
in the Appendix. 



C. Constrained-D calculations 

In Sec. |lT] we have shown that a pathological spill-out 
regime can be triggered by the ferroelectric displacement 
D of the film, as the band offset generally strongly de- 
pends on D. It is therefore important, in order to per- 
form the analysis described in the previous sections, to 
calculate the electronic and structural ground state of a 
metal/ferroelectric interface at different values of D. To 
this end, we can use two different approaches in first- 
principles calculations. The first, and more "traditional" 
approach, involves the construction of capacitor of vary- 
ing thicknesses t, and the relaxation of the corresponding 
ferroelectric ground states within short-circuit boundary 
conditions. Due to the interface-related depolarizing ef- 
fects mentioned in Sec. [ll] (these are strongest in thinner 
films and tend to reduce P from the bulk value Pg), the 
polarization will increase from P = (for t < tent , where 
tcrit is the "critical thickness"^^^jS3.) to P '-^ Ps, in the limit 
of very large thicknesses. This might be cumbersome in 
practice: thicker capacitor heterostructures imply a sub- 
stantial computational cost, due to the larger size of the 
system; this severely limits the range of P values that 
can be studied within short-circuit boundary conditions. 

An alternative, more efficient methodology to explore 
the electrical properties of the interface as a function of 
polarization, is to use the recently developed techniques 
to constrain the macroscopic electric displacement to a 
fixed valuei ^"'^^ With this method, one is able, in princi- 
ple, to access at the same computational cost the struc- 
tural and electronic polarization of the capacitor for an 
arbitrary polarization state. In the specific context of the 
present work, however, there are two drawbacks related 
to the use of the constrained- 1? method as implemented 
in Refs. |33 and [l^. First, fixed- Z? strategies make use 
of applied electric fields to control the polarization of 
the system. This is a problem here, where the metallic- 
ity associated with the space charge which populates the 
ferroelectric film makes such a solution problematic. (If 
a capacitor becomes metallic, it is a conductor and no 
metastable polarized state can be defined at any given 
bias.) Second, our philosophy in this work is to adopt 
"standard" computational techniques, i.e. those that are 
in principle available in any standard electronic structure 
package. 

To this end, we introduce here an alternative way 
of performing constrained- Z? calculations for a metal- 
insulator interface, which does not rely on the direct ap- 
plication of macroscopic electric fields or on the calcu- 
lation of the macroscopic Berry-phase polarization. We 
adopt a vacuum/ferroelectric/metal geometry. To induce 
a given value of the polarization in the ferroelectric film, 
we introduce a layer of bound charges (Q per surface 
unit cell S) at its free surface. If we do so in such a 
way that the surface region remains locally insulating, 
at electrostatic equilibrium, the difference in the macro- 
scopic displacement D on the left and on the right side 
of the surface will exactly correspond to the additional 



13 



surface charge density Q/S. By applying a dipole cor- 
rection in the vacuum region, we ensure that D — Q in 
the region near the surface on the vacuum side; then on 
the insulator side we have exactly 

^-f- (34) 

In practice, the additional charge density is introduced 
by substituting a cation at the ferroelectric surface by a 
fictitious cation of different formal valence. As we are 
interested in exploring intermediate values of I?, we use 
the virtual crystal approximation to effectively induce a 
fractional nuclear charge. 

The reader might have noted that this method to con- 
trol D is just a generalization of Eq. (|T3| to consider 
other forms of "external" charge that are not "free" in 
nature. Indeed, in the most general case, one can state 

V-D(r) =pext(r), (35) 

where D encompasses all bound-charge effects that can 
be referred to the properties of a periodically repeated 
primitive bulk unit, and pcxt contains all the rest (e.g. 
delta-doping layers, metallic free charges, charged ad- 
sorbates, variations in the local stoichiometry, etc.). 
In Eq. (IM)) we simply applied Eq. ([55]) to the vac- 
uum/ferroelectric interface, where the "bound" nature of 
the external charge allows us to control it as an external 
parameter. 

D. Computational parameters 

To demonstrate the generality of our arguments, which 
are largely independent of the fine details of the calcula- 
tion (except for the choice of the density functional) , we 
use two different DFT-based electronic structure codes, 
Lautrec and Siesta42, In both cases, the interfaces 
were simulated by using a supercell approximation with 
periodic boundary conditions4S. A (1 x 1) periodicity of 
the supercell perpendicular to the interface is assumed. 
This inhibits the appearance of ferroelectric domains 
and/or tiltings and rotations of the O octahedra. A ref- 
erence ionic configuration was defined by piling up m 
unit cells of the perovskite oxide (PbTiOa, BaTiOa, or 
KNbOa), and n unit cells of the metal electrode (either a 
conductive oxide, SrRuOa, or a transition metal, Pt). In 
order to simulate the effect of the mechanical boundary 
conditions due to the strain imposed by the substrate, 
the in-plane lattice constant was fixed to the theoretical 
equilibrium lattice constant of bulk SrTiOa (ao = 3.85 A 
for Lautrec and ao = 3.874 A for Siesta). 

To simulate the capacitors in an unpolarized configu- 
ration in Sec. IIV[ we imposed a mirror symmetry plane 
at the central BO2 layer, where B stands for Ti or Nb, 
and relaxed the resulting tetragonal supercells within 
PA/mmm symmetry. For the ferroelectric capacitors de- 
scribed in Sec. [via- second minimization was carried out. 



with the constraint of the mirror symmetry plane lifted. 
Tolerances for the forces and stresses are 0.01 eV/A and 
0.0001 eV/A^, respectively. Other computational param- 
eters, specific to each code, are summarized below. 



1. Lautrec 

Calculations in Sec. IIV Bl and IV Al wcre performed with 
Lautrec, an "in- house" plane- wave code based on the 
projector-augmented wave method4^ We used a plane- 
wave cutoff of 40 Ry and a 6 x 6 x 1 Monkhorst-Pack^Si^ 
mesh. As the systems considered here are metallic, we 
adopted a Gaussian smearing of 0.15 eV to perform the 
Brillouin-zone integrations. 



2. Siesta 

Computations in Sec. IIV Al and IVBI on short-circuited 
SrRuOa/PbTiOs and SrRuOs/BaTiOa capacitors were 
performed within a numerical atomic orbital method, as 
implemented in the Siesta code."'^ Core electrons were 
replaced by fully-separable^ norm-conserving pseudopo- 
tentials, generated following the recipe given by Troullier 
and Martinsj^ Further details on the pseudopotentials 
and basis sets can be found in Ref. HI- 

A 6 X 6 X 1 Monkhorst-Pack^Si^ mesh was used for the 
sampling of the reciprocal space. A Fermi-Dirac distri- 
bution was chosen for the occupation of the one-particle 
Kohn-Sham electronic eigenstates, with a smearing tem- 
perature of 0.075 eV (870 K). The electronic density, 
Hartree, and exchange-correlation potentials, as well as 
the corresponding matrix elements between the basis or- 
bitals, were computed on a uniform real space grid, with 
an equivalent plane-wave cutoff of 400 Ry in the repre- 
sentation of the charge density. 



IV. RESULTS: PARAELECTRIC CAPACITORS 
A. Non-pathological cases 

In the centrosymmetric unpolarized reference struc- 
ture, some metal/ferroelectric interfaces such as 
BaTiOa/SrRuOa or PbTiOa/SrRuOa are "well-behaved" 
within LDA. [We focus here on the Ti02/SrO termina- 
tion - the properties of the alternative (Ba,Pb)0/Ru02 
termination might differ.] This conclusion emerges from 
the analysis shown in Fig. [6]for the PbTiOa-based capac- 
itor; qualitatively similar results, not shown here, are ob- 
tained for the BaTiOa-based capacitor. Figure [B^a) rep- 
resents schematically the Schottky barriers for electrons 
{(fin) and holes {(j)p) at the ferroelectric/metal interfaces, 
computed using the nanosmoothed electrostatic poten- 
tial method described in Sec. Ill Al The bottom of the 
conduction band of the ferroelectric lies above the Fermi 
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FIG. 6: (a) Schematic representation of and (j>p in a-n un- 
polarized SrRuOs/PbTiOs/SrRuOs capacitor. Ey, Ec, Ef 
and A(V^) were defined in Sec. Ill Al The calculated values are 
also indicated in the Figure. The black solid line represents 
—Vu{z). The dashed line represents the hypothetical position 
of the CBM if Ec were shifted to reproduce the experimen- 
tal band gap. (b) Profile of pfrco as defined in Eq. (|26|l . (c) 
Profile of the layer- by-layer polarization Pf . The size of the 
capacitor corresponds to n = 5.5 unit cells of SrRuOa and 
m = 12.5 unit cells of PbTiOa. Only the top half of the 
symmetric supercell is shown. 



level of the metal (0„ amounts to 0.38 eV for the PbTiOa- 
based capacitor, and only to 0.19 eV in the BaTiOa-based 
case). Note that, if the experimental band gap could be 
reproduced in our simulations, would be much larger 
[dashed lines in Fig. [BJa); we have taken the experimen- 
tal indirect gap of the cubic phase of PbTiOs, 3.40 eV^^ 
and assumed that the quasiparticle correction on the va- 
lence band edge is negligible] . The results summarized in 
Table HIl indicate that, in all the cases discussed here, dif- 
ferent methodologies yield Schottky barrier values that 
are consistent within a few hundredths of an eV. The 
flatness of the profile of the nanosmoothed electrostatic 
potential at the central layers of PbTiOs confirms the ab- 
sence of any macroscopic electric field, as expected from 
a locally charge-neutral and centrosymmetric system. 

Figure ini^b) displays pheciz), as defined in Sec. IIII B"T] 
As expected, phcdz) has a rapid decay in the insulat- 
ing layer, consistent with the evanescent character of the 
metallic states (MIGS): these cannot propagate in the 
insulator as their energy eigenvalues fall within the for- 
bidden band gap. Fig. [SJc) shows the layer-by-layer po- 
larization, , computed using Eqs. Consis- 
tent with the absence of space charge, the profile is 
remarkably flat. Due to the imposed mirror-symmetry 
constraint, P^ also vanishes inside the ferroelectric ma- 
terial. 

Fig. [7] shows the layer-resolved PDOS of the Ti(3s) 
semicore peaks, the 0(2s) peak, the upper valence band 
and the lower conduction band (black curves, shaded 
in gray). On top of the heterostructure PDOS we su- 



TABLE IL LDA values of (j)„ and (jjp, obtained with two 
different methods: using the decomposition into Ey^c and 
A{V) (BS -I- Lineup), or using the method of Sec. IIII A 31 
(Semicore). In the "Semicore" case we used the sharp Ti(3s)- 
derived peak of the LDOS (extracted from the central Ti02 
layer of the capacitor) to align the energies of the bulk band 
edges with the supercell Fermi level. 



Capacitor 



BS + Lineup Semicore 



SrRuOa/PbTiOs/SrRuOs 

(f)p (eV) 0.97 0.99 

(j>„ (eV) 0.38 0.37 

SrRuOs/BaTiOs /SrRuOa 

(j}p (eV) 1.39 1.41 
4>„ (eV) 019 0.23 




-18 -16 
Energy (eV) 

FIG. 7: (Color online) PDOS of the inequivalent Ti02 layers 
in the unpolarized PbTiOa/SrRuOa capacitor (solid curves 
with gray shading). The bottom curve lies next to the elec- 
trode, the top one lies in the center of the PbTiOs film. Only 
the PDOS on half of the symmetric supercell are shown. The 
bulk PDOS curves (red dashed) are aligned to match the 
Ti(3s) peak at iJ ~ —57 eV. The Fermi level is located at 
zero energy. 



perimpose the bulk PDOS, calculated with an equiva- 
lent /c-point sampling and aligned with the Ti(3s) peak 
(dashed red curves). Note that all PDOS curves were 
calculated using Eq. (22), and the smearing function gpD 
of Eq. (|23bp with a = 0.075 eV, consistent with the 
parameters used in the calculation. The PDOS of the 
conduction and valence bands converges fairly quickly 
to the bulk curve when moving away from the interface 
- they are practically indistinguishable already at the 
fourth layer. The estimated energy locations of the con- 
duction and valence bands converge even faster [these are 
directly related to the shifts of the Ti(3s) state, which are 
less affected by quantum confinement effects]. All curves 
except those adjacent to the electrode interface vanish at 
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FIG. 8: (Color online) LDOS integrated over the Nb02 lay- 
ers of the KNbOa/SrRuOa heterostructure (solid curves with 
gray shade). The bottom curve lies next to the electrode, 
the top one lies in the middle of the KNbOa film. Only the 
PDOS on half of the symmetric supercell are shown. The bulk 
LDOS (red dashed curves) are aligned to match the valence 
and conduction band edges. The Fermi level is located at zero 
energy. 

the Fermi level, confirming the absence of charge spill-out 
in this system. 

As a summary of this Section we can conclude that, 
when a centrosymmetric unpolarized interface is non- 
pathological in the sense that the bottom of the conduc- 
tion band of the ferroelectric is above the Fermi level of 
the metal, (i) the free charge, as defined in Sec. HUB 1[ 
vanishes due to the absence of charge spill-out; (ii) the 
local polarization profile (Sec. IIII B 31) is perfectly fiat 
as the interface-induced polar lattice distortions heal 
rapidly (within the first unit cell); (ni) the LDOS/PDOS 
vanishes at the Fermi level, except for one or two inter- 
face layers where the signatures of the MIGS might be 
still present (they are barely detectable in the curves of 
Fig. El). 



B. Pathological cases 

We analyze now two examples of capacitors that 
are characterized by a pathological band alignment 
already in their centrosymmetric reference struc- 
ture: Nb02-tcrminatcd KNbOs/SrRuOs, and Ti02- 
terminated BaTiOg/Pt. This choice of materials is mo- 
tivated by the fact that there exist recent theoretical 
works on these systems^^i^ where the consequences of 
the pathological band alignment were neglected. 



1. KNbOs/SrRuOs 

We construct a heterostructure consisting of m=6.5 
KNbOa unit cells and n=7.5 SrRuOa cells, for a total 




FIG. 9: (Color online) Calculated free charge for paraelec- 
tric SrRuOs/KNbOa/SrRuOs heterostructure. Black curve: 
planar-averaged pfjpc- Red dashed: Pfroc i ns^riosmoothed using 
a Gaussian filter. Blue symbols: finite differences of the local 
Pj (shown as a black curve in Fig. llOp . calculated using the 
Wannier-based layer polarization described in Sec. IIII B 21 



of 14 perovskite units; we use symmetrical Nb02 (SrO) 
terminations of the KNbOa (SrRuOa) film. After full 
relaxation with a mirror symmetry constraint at the cen- 
tral Nb02 layer, we perform the analysis of the LDOS, 
the conduction charge and the local polarization as ex- 
plained in Sec. IIIII In Fig. |8] we show the local density of 
states integrated over the Nb02 layers (the bottom one 
is adjacent to the electrode interface, the top one lies on 
the mirror plane in the middle of the film) . The unphys- 
ical Ohmic band alignment is evident from the location 
of the conduction band bottom ~ the whole film is clearly 
metallic. This points to the pathological situation that 
is sketched in Fig. [3l Note that the LDOS does not con- 
verge to the bulk curve anywhere in the heterostructure. 
There are non-trivial shifts of all peaks that make it dif- 
ficult to identify a well-defined alignment with the bulk 
curves. In Fig. [8] we choose to align the 0(2s)-derived 
feature at i? ~ — 19 eV. In this specific system, align- 
ing the 0(2s) peaks appears to yield a reasonably good 
match of the conduction and valence band edges (the 
most relevant features from a physical point of view); 
this, however, leads to a marked mismatch, e.g. in the 
position of the semicore Nb(4s) state. We show in the fol- 
lowing that these effects stem from a number of (rather 
dramatic) electrostatic and structural perturbations act- 
ing on the KNbOa film, which are a direct consequence 
of the pathological band alignment. 

First we show that the non-vanishing LDOS at the 
Fermi level results in a sizeable spill-out of conduction 
charge into the ferroelectric film. To that end, we plot 
Pfj.(,p(z), which represents the planar average of the arti- 
ficially populated part of the KNbOa conduction band, 
and the corresponding nanosmoothed version, 'pfj.^^(z), 
in Fig. [3] respectively as black continuous and red dashed 
lines. The additional electron density in the ferroelec- 
tric region is apparent, and reaches a maximum of about 
0.15 electrons in the central perovskite unit cell. Such 
a density is significant - it can be thought as resulting 
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FIG. 10: (Color online) Local polarization profile in the 
SrRuOa/KNbOa/SrRuOs capacitor. Black circles: polariza- 
tion from Wannier-based layer polarizations. Red squares: 
approximate polarization from "renormalized" Born effective 
charges (see Sec. IIIIB 3|) . Analogous results for a paraelectric 
SrRuOa/BaTiOs/SrRuOs capacitor are shown for compari- 
son (blue diamonds). 

from an unrealistically large doping of, e.g. one Sr^+ 
cation every six or seven ions. However, unlike in a 
doped perovskite, the spurious electron spill-out here is 
not compensated by an appropriate density of heterova- 
lent cations. The system is therefore not locally charge 
neutral, and as a consequence strong, non-uniform elec- 
tric fields arise in the insulating film that act on the ionic 
lattice. 

In order to elucidate how the underlying polarizable 
material responds to such an electrostatic perturbation, 
we plot in Fig. [10] the effective polarization profile in 
the KNbOs film calculated in two ways, (i) the rigor- 
ous Wannier-function analysis of the layer polarizations 
and (ii) the approximate expression based on the renor- 
malized bulk dynamical charges. The matching between 
the curves is excellent, indicating that the approximate 
Z*-based formula provides a reliable estimate of P{z); 
this suggests that the electrostatic screening is indeed 
dominated by structural relaxations, as anticipated in 
Sec. ini and as expected in a ferroelectric material. To 
substantiate this point, we compare in Fig. [TT] the re- 
laxed layer rumplings in KNbOs/SrRuOa to those of 
the non-pathological case, PbTiOa/SrRuOa, discussed in 
Sec. IIV Al The KNbOa film is characterized by strong 
non-homogeneous distortions, which are consistent with 
the polarization pattern shown in Fig. 1101 Conversely, 
the distortions are negligible in the PbTiOs/SrRuOa ca- 
pacitor, where all the oxide layers are essentially fiat. 

The polarization profile Pj is characterized by a uni- 
form, negative slope. This nicely confirms the prediction 
of our semiclassical analysis in Sec. |TT]of a uniform linear 
decrease of D{z) throughout the film. Pj varies from 0.3 
to -0.3 C/m^ when moving from the bottom to the top in- 
terface. Note that such spatial variation is completely ab- 
sent in, e.g., isostructural paraelectric BaTiOa/SrRuOs 
(diamonds in Fig.UHl), and PbTiOa/SrRuOg [Fig. M^c)] 
capacitors, where the profile is remarkably fiat with P 




FIG. 11: (Color online) Layer rumplings (cation-oxygen ver- 
tical relaxations) in the centrosymmetric KNbOa/SrRuOs 
(black line, empty circles) and PbTiOa/SrRuOa (red line, 
filled circles) capacitors. Dashed vertical lines indicate the 
location of the BO2 planes. The shaded areas correspond to 
the SrRuOs electrode region. 

vanishing throughout the film. We stress that the non- 
uniform perturbation experienced by KNbOs/SrRuOa 
is qualitatively different from a ferroelectric distortion, 
which involves an almost perfectly rigid displacement of 
the ionic sublattices: in absence of space-charge effects, 
a macroscopically uniform rumpling pattern across the 
film is typically found 1^ 

To demonstrate that the spatial variation in P{z) is 
directly related to phee according to Eq. (IT^ . we per- 
form a numerical differentiation of the polarization pro- 
file derived from the Wannier-based layer polarizations. 
The result, plotted in Fig. [S] as a blue line, shows an 
essentially perfect match between dP/dz and — pftoc il- 
lustrating the fact that the polarization profile is really a 
consequence of KNbOa responding to the spurious pop- 
ulation of the conduction band, rather than of interface 
bonding effectsi^ 



2. BaTiOa/Pt 

We next present results of an analogous investiga- 
tion for a paraelectric (BaTi03)„i/(Pt)„ capacitor, with 
m — 8.5 and n ~ 11. We consider symmetric Ti02 
terminations, with the interfacial O atoms in the on- 
top positions. (Note that this interface structure is dif- 
ferent than the AO-terminated films simulated, e.g. in 
Refs. [13 and [H, where a Schottky-like band offset was 
found.) We find this interface to have a pathological band 
alignment, similar to the KNbOa/SrRuOa case discussed 
above. The comparative analysis of the bound-charge 
polarization profile and of the excess conduction charge, 
shown in Fig. I12[ again shows excellent agreement be- 
tween pboc{z) and the compensating bound charge. The 
effect is analogous to KNbOa/SrRuOa, with an overall 
magnitude which is smaller by roughly a factor of two; 
the polarizations at the two extremes of the film reach 
values of about ±0.15 C/m^. 

The almost perfect similarity in behavior between 
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FIG. 12: (Color online) (a) Calculated free charge and (b) 
local polarization profile for a paraelectric Pt/BaTiOa/Pt ca- 
pacitor with Ti02-type interfaces. All symbols have the same 
meaning as in Fig. [U] and Fig. 1101 

these two chemically dissimilar systems is further proof 
that the unusual effects described here and in Ref. [sl 
- the apparent head-to-head domain wall in the ferro- 
electric film - have little to do with the bonding at the 
interface, but are merely a consequence of the artificial 
charge spill out, as discussed in Sec. [TTl 

Before moving on to the next Section we briefiy com- 
ment on the physical nature of the conduction charge 
that spills into the ferroelectric film. In particular, it is 
important to clarify that the charge densities plotted in 
Fig. [S] and Fig. [T^ a) indeed originate from population 
of the conduction band of the insulator, and not from 
metal-induced gap states (MIGS) as some authors have 
recently argued»^ First, all charge density plots show a 
maximum in the middle of the ferroelectric layer, rather 
than a minimum, which one would expect if the former 
hypothesis were true, given the evanescent character of 
the MIGS. Second, if MIGS were present they would be 
clearly identifiable in the local density of states; however, 
the LDOS plotted in Fig. [5] shows no evidence of quan- 
tum states lying within the energy gap of the KNbOa 
film. Therefore, we must conclude that these are gen- 
uine conduction band states, and not MIGS. The maxi- 
mum of pfrcc in the middle of the ferroelectric film can be 
interpreted either as a quantum confinement effect [the 
lowest-energy solution of the electron-in-a-box problem 
is indeed a sine function with a shape reminiscent of the 
Pfree plots of Fig. [3] and Fig. fT^ a)]. and/or as a result of 
the dielectric nonlinearities discussed in Sec. Ill C 31 

C. Estimating the "pre-spill" band offset 

We mentioned in Sec. |lT] that, whenever an elec- 
trode/ferroelectric interface enters the pathological spill- 
out regime, the transfer of charge into the conduction 
band of the insulator produces an upward shift of the 
CBM. This effect prevents a direct, unambiguous de- 
termination of the interface parameter 0^. To circum- 
vent this problem, and obtain an approximate estimate 
of the negative "pre-spill" Schottky barrier 0°, we use 



FIG. 13: (Color online) n-type Schottky barrier as a func- 
tion of interface doping in KNbOa /AO-terminated SrRuOs, 
where A is a fictitious atom with atomic number Z — 37 + x. 
Only the Sr atoms at the interfacial layer are replaced by this 
fictitious atom. The dashed line is a linear regression of the 
data between a; = and x — 0.3, where the interface is non- 
pathological from the band alignment point of view. Blue 
and red empty symbols represent, respectively, the results for 
X = 0.5 and x = 1.0, where the interface is already patholog- 
ical. All values were obtained from Eq. H24p . using either the 
Nb(4s) (squares) or the 0(2s) (circles) semicore peaks of the 
central Nb02 layer as a reference. 
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FIG. 14: (Color online) (a) Conduction electrons, and (b) 
local (Wannier-based) polarization profiles extracted from the 
calculations with x =0.0, 0.1, 0.2 and 0.3 (filled circles, thin 
black curves); 0.5 (empty blue circles, dashed blue curve) and 
1.0 (empty red circles, solid red curve). In (a) only half of the 
KNbOa film is shown. 

an approach inspired by a recent worki^ The authors of 
Ref. [53 show that the Schottky barrier at the interface be- 
tween a perovskite insulator (SrTiOa) and a perovskite 
electrode (Lao.yAo.aMnOa, where A is Ca, Sr, or Ba) 
evolves linearly as a function of the compositional charge 
of the interface layer. (This interface layer is of the type 
Laa;Sri_2;0, where x interpolates between a -1-3 and a +2 
cation.) Of course, this linear behavior refers to a range 
of X values where the interface is non-pathological; our 
arguments indicate that as soon as the system enters the 
spill-out regime, the value of saturates to a nearly 
constant value. Based on this observation, if one knows 
the linear behavior of (j)n in a range of x values for which 
the interface is non-pathological, one can extrapolate this 
straight line to the values of x which cannot be directly 
calculated, and obtain an estimate for 0^. 

We apply this strategy to the same KNbOs/SrRuOs 
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capacitor system described in Sec. IIVB II To tune the 
interface charge, we replace the Sr cation in the inter- 
face SrO layer with a fictitious atom of fractional atomic 
number Z = 37 + x — 1 corresponds to the example 
already shown in Sec. IIVB 11 with a charge-neutral SrO 
interface layer, and x = corresponds to a RbO layer of 
net formal charge -1. The results for the Schottky barrier 
are plotted in Fig.[T31 The region from x = 0.0 to x = 0.3 
is non-pathological and shows an almost perfectly linear 
evolution of (dashed line) . By extrapolating this lin- 
ear trend to a; = 1 we obtain (/)„ ~ —1.2 eV, which is 
about 1 eV lower than the value calculated from first 
principles. This confirms the remarkable impact of the 
space-charge effects described in Sec. IIIC II Assuming a 
polarization of ^ 0.3 C/m^ for KNbOs near the inter- 
face, a potential drop of 1 eV would be accounted for 
by an effective screening length of 0.3 A at the electrode 
interface. This value is quite reasonable, and similar in 
magnitude to those reported in Table U . 

In order to examine the crossover between the Schottky 
(non-pathological) and the Ohmic (pathological) regimes 
in terms of the analysis tools developed in this work, we 
plot in Fig. [13] the polarization profiles and the density 
of conduction electrons for each of the calculations sum- 
marized in Fig. [T31 These plots confirm that from x — 
to X = 0.3 the capacitors are non-pathological, with ab- 
sence of conduction charge in the insulating region [panel 

(a) , thinner lines] and a flat polarization profile [panel 

(b) , filled circles - all these curves overlap on this scale]. 
Conversely, at a; = 0.5 the conduction band starts popu- 
lating significantly [thicker dashed blue line in (a), empty 
blue circles in (b), note that the corresponding points 
in Fig. [T2] starts to depart from the linear regime]. At 
a; = 1.0 population of the conduction band has become 
dramatic, and so is the corresponding slope in the polar- 
ization profile. The departure from linearity in Fig. [T3] 
is correspondingly large. Note that the use of either the 
Nb(4s) or the 0(2s) semicore peaks in Eq. ([M)) yields 
identical results in the non-pathological regime (the filled 
squares and circles overlap in Fig. [T51) . Conversely, the 
result depends significantly on this (completely arbitrary) 
choice at X = 0.5, and even more so at x = 1.0 (the cir- 
cles and squares split). This is another proof that in the 
pathological regime the band lineup is ill-defined - due to 
the electrostatic effects discussed throughout this work, 
the LDOS does not converge to a bulk-like value in the 
center of the KNbOa film (see Fig. [8]), and there is no 
obvious reference energy to determine the offset. 



V. RESULTS: FERROELECTRIC CAPACITORS 

As discussed in the Introduction, although some 
of the unpolarized reference structures (e.g. the 
PbTiOs/SrRuOs interface) appear artifact-free within 
LDA, because of the strong dependence of the Schot- 
tky barrier on D [Eq. [5] they might become problematic 
when the constraint of mirror symmetry is lifted and the 
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FIG. 15; (Color online) Results for the polarized 
PbTiOa/SrRuOa interface for increasing polarization of the 
film, (a) planar averaged Pfrec- Black, red, green and blue 
curves correspond to the results for d =0.20, 0.40, 0.60 and 
0.74 e, respectively. The sharp peaks in Pf^ee correspond to 
the Ti ions in the PbTiOs film, (b) layer polarizations from 
the Wannier-based analysis. Same color code as in (a). 



system is polarized. To address this issue, in this section 
we first use the fixed-D strategy described in Sec. IIII CI 
to explore the behavior of the PbTiOs/metal interface 
over a wide range of polarization states. Then, we will 
demonstrate that the behavior that we calculate using 
the fixed- method corresponds exactly to that of a true 
short-circuited capacitor by performing more "standard" 
large-scale calculations for a few selected thickness val- 
ues. 



A. Open-circuit calculations 

We construct a vacuum/PbTiOa/SrRuOa heterostruc- 
ture as explained in Sec. IIII CI The reduced macroscopic 
displacement field^-, d = DS, is controlled by substi- 
tuting the Ti at the PbTiOs/vacuum interface with a 
fictitious cation of atomic number Zi^ft — AQ -\- d (i.e. 
Zr for d=0). The thickness of the PbTiOs slab is set to 
5 unit cells, and that of SrRuOa to 4; other computa- 
tional parameters are kept the same as in the rest of this 
work. We considered four different values of d: 0.2, 0.4, 
0.6 and 0.74, the latter one corresponding to the ferro- 
electric ground state of PbTiOs at the SrTiOa in-plane 
lattice constant. In each case, we verify by examining 
the LDOS that the free surface remains locally insulat- 
ing; therefore, the macroscopic D ~ d/S in the film cor- 
responds exactly to the value enforced by the artificial 
pseudopotential. 

The evolution of p^cc and of the Wannier-based layer 
polarization profiles for 0.2 < d < 0.74 is shown in 
Fig. [13 It is apparent from the plots of pfrec that already 
for the smallest value of the polarization [d = 0.20, black 
curve in Fig.fTSTa')] the Ti02 layer closest to the electrode 
has an important density of conduction electrons. This 
is expected, as the evanescent tails of the metal-induced 
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FIG. 16: (Color online) Calculated results for the fully polar- 
ized PbTiOa/SrRuOa interface at d = 0.74. (a) local polar- 
ization from the Wannier-based layer polarizations, and (b) 
planar averaged pf^gj. (black curve) , macroscopically averaged 
PfroG (■'S'i dashed curve), and finite differences of the polariza- 
tion shown in the panel (a) (blue squares) for a,n m — 8-unit 
cell thick PbTiOa film. Panels (c) and (d) are the correspond- 
ing figures for an ?n = 5-unit cell thick PbTiOa film. The 
sharp peaks in Ptroc correspond to the Ti ions in the PbTiOs 
film. 



gap states (MIGS) penetrate into the insulating region 
for some distance at any metal/insulator junction. How- 
ever, these states do not propagate very far, and already 
at the second Ti02 layer they are barely noticeable on 
the scale of the plot. At d = 0.4 [red curve in Fig. [TSl a)] 
the peak on the second Ti02 layer significantly increases 
in magnitude, and a new small peak appears at the third 
Ti02 layer. Analysis of the local density of states (not 
shown) shows that these new peaks are conduction band 
states of PbTiOs, rather than evanescent SrRuOa states. 
The reason why ptrce decays relatively fast when mov- 
ing into the insulator is due here to the internal field in 
PbTiOs, which generates a confining wedge-like poten- 
tial. We stress again that this mechanism is fundamen- 
tally different from the usual quantum-mechanical damp- 
ing of the MIGS that fall in a forbidden energy window of 
the insulator. We identify this mechanism with the on- 
set of the Schottky breakdown, which becomes increas- 
ingly apparent if the polarization of the film is further 
increased to d — 0.60 [green curve in Fig. [TSfa).] As 
in the discussion of the paraelectric capacitors, the pres- 
ence of the space charge is reflected in the progressive 
"bending" of the layer polarization profile [Fig. [TSlb)]. 

At d = 0.74, the population of the conduction band be- 
comes rather dramatic, and the charge distributes over 
the whole film. Here, the space charge is no longer con- 
fined by the depolarizing field: in the fully polarized 
ferroelectric state the internal field of PbTiOs vanishes. 
Therefore, the intrinsic carriers are only loosely bound 



to the interface by the band bending effect, analogous to 
the mechanism that confines the compensating carriers 
at the LaAlOs/SrTiOs interfacei^ Since the dielectric 
permittivity of PbTiOs is rather large, the band bend- 
ing is very efficiently screened, and the distribution of 
charge can reach quite far into the insulator. To demon- 
strate this fact, we have repeated the simulation with the 
same value of d—Q.7A, but with a thicker PbTiOs film of 
8 unit cells [Fig. [16] (b)]; indeed, the conduction electrons 
redistribute over the whole volume of the film to mini- 
mize their kinetic energy. Thus, in close analogy to the 
LaAlOs/SrTiOs case, the metallization of the fully po- 
larized PbTiOs film at (i=0.74 can be thought as a form 
of "electrostatic doping" induced by spill-out of electrons 
from the electrode to the PbTiOs conduction band. We 
shall further elaborate on this point in Sec. IVI Fl 

Fig- llSf b) illustrates a further important consequence 
of the charge spill-out regime, which was mentioned al- 
ready in Sec. lII C 2l in the pathological regime the dipoles 
that lie closest to the electrode interface may appear 
"pinned" to a fixed value. This is indeed the case for the 
Ti02 layer adjacent to the electrode, which seems to sat- 
urate at ~ 0.08 nC/m for increasing values of D. Again, 
we caution against interpreting this dipole pinning effect 
as a robust physical result. 



B. Short-circuit calculations 

To demonstrate in practice that the conclusions of 
Sec. lV Al inferred by using open-circuit boundary condi- 
tions, are directly relevant to short-circuited capacitors, 
we have performed simulations on SrRuOs /PbTiOs het- 
erostructures, with m = 12.5 and n = 5.5. A soft-mode 
distortion of the bulk tetragonal phase, inducing a polar- 
ization perpendicular to the interface, is superimposed on 
the PbTiOs layers of the previous unpolarized configura- 
tions discussed in Sec. IIV Al Then the atomic positions 
of all the ions, both in the ferroelectric and in the metallic 
electrodes, and the out-of-plane stress are relaxed again 
with the same convergence criteria as before. 

By means of the approximate Eq. p3p . derived in Sec. 
nil B 3[ we computed the local layer-by-layer polarization, 
Pf, plotted in Fig. [IT] (a). Far enough from the inter- 
face, the polarization profile is rather uniform, with a 
polarization that amounts to 0.53 C/m^ in PbTiOs (64 
% of the strained bulk polarization) , which we identify as 
the macroscopic P of the PbTiOs film. This corresponds 
to d ~ 0.5, i.e. a value that in the open-circuit study of 
the previous section (see Fig. [T5|) we found to be already 
pathological. To verify that the same happens here, we 
analyze the density of conduction electrons. The planar 
average of ptrcc(r) for the relaxed polar configuration is 
plotted in Fig. [TTTb). The existence of a charge populat- 
ing the Ti 3c? orbitals is evident from the peaks of pfroc(-z) 
at the Ti02 layers, which are detectable up to four unit 
cells away from the interface. Indeed the profile of the 
conduction charge appears to be intermediate between 
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FIG. 17: (Color online) (a) Profile of the layer- by- layer po- 
larization Pf , defined in Eq. (|33|) . in the relaxed polar con- 
figuration of a short-circuited SrRuOa/PbTiOa/SrRuOa ca- 
pacitor. The dashed line represents the bulk spontaneous po- 
larization under the same in-plane strain as in the capacitor, 
(b) phee{z) as defined in Eq. (|26[) (black solid line), and its 
nanosmoothed average pfrcc(z) (red dashed line). The blue 
line represents the profile of the bound charge, computed as 
a finite-difference derivative of Pf . 



the d = 0.4 and d = 0.6 cases of Sec. IV Al consistent 
with the present estimate d ^ 0.5. 

As we already anticipated in the previous Sections, 
Pfrco is responsible for non-trivial lattice relaxations, 
which act to screen the electrostatic perturbation. Fig. 
[TTT a) indeed shows a small bending of the local polariza- 
tion profile, starting roughly three unit cells away from 
the top interface and with a negative slope of the local po- 
larization, . To prove that such a spatial variation of 
P{z) [which in PbTiOs provides a reasonably accurate es- 
timate of the local electric displacement D{z)] is a direct 
consequence of the presence of the non- vanishing conduc- 
tion charge [recall Eq. (jl3p ]. we numerically differentiati- 
ate the polarization profile and compare it with p{,-(.c{z) 
in Fig. [T7i;b). As in the SrRuOa/KNbOa/SrRuOa unpo- 
larized case (see Fig. 19]) , the bound charge (divergence of 
P) accurately neutralizes the conduction charge (nanos- 
moothed profile of Pfrco)- 

In order to further prove that the present case fits into 
the model description of Sec. Ill C 21 in Fig. [18] we plot 
the layer-resolved PDOS. The curves were constructed 
exactly as in Fig. [71 except that (i) the capacitor is 
now polarized; and (ii) consistent with the discussion of 
Sec. nil A 31 we set up the bulk reference calculation by 
using the PbTiOs structure extracted from the polarized 
supercell calculation (i.e. with atomic distortions and 
out-of-plane strain consistent with a polarization of 0.53 
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FIG. 18: (Color online) Layer-by-layer PDOS on the Ti02 
layers of the polar SrRuOs/PbTiOa/SrRuOs ferroelectric ca- 
pacitor. Meaning of the lines as in Fig. [T] but now the PDOS 
on all the Ti02 layers are plotted (there is no longer a mirror 
symmetry plane). The squares represent the position of the 
local band edges, computed following the recipe of Sec. lIIIiO] 
The dashed lines are a linear interpolation of the calculated 
local band edges. 



C/m^). The agreement is again very good, showing that 
our approximation of neglecting the macroscopic depolar- 
izing field in the bulk reference calculation is a reasonable 
one, and that the most important effects on the LDOS 
originate from the lattice distortions. In the capacitor 
we clearly distinguish two regions. In the lower part of 
the PbTiOs film, the PDOS at the Fermi level vanishes, 
which implies that the system is locally insulating. Fur- 
thermore, the PDOS in each layer appears rigidly shifted 
with respect to the neighboring two layers, consistent 
with the presence of a depolarizing field. In the upper 
region, close to the top electrode, the PDOS crosses the 
Fermi level and the system is locally metallic. All these 
features are in full agreement with the scheme drawn in 
Fig.m 

In Fig. [TSj we also plot the estimated band edges for 
each layer, Ey^cU), obtained from Eq. We used 

the semicore Ti(3s) peak at each layer as Esdj), and 
we calculated the bulk contributions in Eq. ([Ml) from a 
non-self-consistent bulk calculation (based on the ground 
state charge density of the bulk reference calculation at 
P =0.53 C/m^ described above) that included the high- 
symmetry fc-points. The resulting data points lie very ac- 
curately on a straight line. By extrapolating this straight 
line, we see that it crosses the Fermi level near the fourth 
PbTiOs cell from the top electrode interface. This illus- 
trates the pathological character of the band alignment 
in this system, consistent with the model of Fig. SI 

As a final remark, we mention that we performed sim- 
ilar calculations for a polarized BaTiOa/SrRuOs capaci- 
tor, and found a very similar scenario, with the conduc- 
tion band locally crossing the Fermi level as the mirror 
symmetry plane is lifted and a spontaneous polarization 
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is allowed to develop. In general, the onset of such a 
pathological regime has important consequences on many 
physical properties of the capacitor, as we shall discuss 
in the following Section. 

VI. DISCUSSION 

In this Section we discuss the important aspects of our 
work in the context of the existing literature. The dis- 
cussion is organized in several categories, corresponding 
to the different properties of a ferroelectric/electrode in- 
terface (or, more generally, of a perovskite material) that 
might be affected by the (more or less spurious) presence 
of free charges in the system. 




FIG. 19: (Color online) Schematic representation of the im- 
pact of the charge spillage on the ferroelectric stability. M 
are the metal electrodes, and FE is the ferroelectric film. The 
polarization points to the right. 



A. Structural properties of the film 

The authors of Ref. [H studied KNbOa thin films 
placed between symmetric metallic electrodes (either 
SrRuOa or Pt) under short-circuit electrical boundary 
conditions. In the SrRuOa case, the layer-by-layer po- 
larization pointed in opposite directions at the top and 
bottom interfaces for all thicknesses, creating 180° head 
to head domain walls, which were denominated interface 
domain walls (IDW). The physical origin of the IDW 
was attributed to a strong bonding between interfacial 
Nb and O atoms, which would induce a "pinning" of the 
interface dipoles to opposite values at the top and bottom 
electrode interfaces. 

Here we have demonstrated with analytical derivations 
and practical examples that both the inhomogeneous po- 
larization and the "dipole pinning" effect are clear sig- 
natures of a pathological band alignment. In particu- 
lar, in an unpolarized KNbOs/SrRuOa capacitor anal- 
ogous to those simulated by Duan et al.^^, we obtain a 
monotonously decreasing polarization profile, from (^0.3 
C/m^) at the bottom interface to an opposite value of 
^-0.3 C/m^ at the top, in excellent agreement with the 
results of Duan and coworkers.- - In contrast with the 
conclusion of Ref. [sgI , however, here we find that the mi- 
croscopic origin of this strong inhomogeneous polariza- 
tion is the spillage of charge from the metallic electrode 
to the bottom of the conduction band of KNbOs, rather 
than a bonding effect. 

These findings have important consequences concern- 
ing the physical understanding of the system with regard 
to the relevant observables. First, the ferroelectric ma- 
terial becomes in fact a metal, and such a device would 
respond Ohmically with a large direct DC current that 
would make switching difficult or impossible. This ques- 
tions the appropriateness of interpreting the "average" 
polarization of the film as a macroscopic physical quan- 
tity that can be measured in an experiment (see next 
Section). Second, our arguments indicate that two es- 
sential factors governing the equilibrium free charge dis- 
tribution (and hence the spatial variation of P) are the 



conduction band structure (e.g. the density of states) 
of the ferroelectric material, and the interface band off- 
set. Both ingredients are absent in traditional Landau- 
Ginzburg models, e.g. those used in Ref. [s^ to inter- 
pret the above data on KNbOa/SrRuOa capacitors, or 
in Ref. [s^ to interpret qualitatively similar results for a 
hole-doped BaTiOa/SrRuOs interface. We therefore cau- 
tion against overinterpreting the results of such models, 
as they might fail at capturing the relevant physics of 
the free-charge equilibration. A promising route towards 
overcoming these limitations appears to be the model 
Hamiltonian approach proposed in Ref. [s^ . Extending 
that strategy to the case of a metal/ferroelectric inter- 
face will be an interesting subject of further research. 



B. Stability of the ferroelectric state 

The pathological spill-out of charge has important con- 
sequences on the spontaneous polarization of a ferroelec- 
tric capacitor. To give a qualitative fiavor of such an 
effect, we consider the case of a capacitor that is only 
partially metallic, i.e. there is a depolarizing field that 
keeps the carriers confined to the pathological side as 
sketched in Fig. [TW a). We further consider two symmet- 
ric electrodes, i.e. characterized by identical values of (/)° 
(that we assume positive) and XeS- Assuming a mon- 
odomain state, there are then two stable configurations, 
related by a mirror symmetry operation. As (j)'^ is posi- 
tive, upon application of an electric field there will always 
be an insulating region in the middle of the film, i.e. the 
polarization can be switched without passing through a 
globally metallic state. 

To appreciate the impact of the charge spill-out on the 
spontaneous polarization of the film, it is useful to look 
at the schematic band diagram of Fig. [TW a). where the 
conduction band bottom goes below the Fermi level in 
proximity of the right electrode (red area). This induces 
metallicity in a significant portion of the film (light grey 
shaded area, up to the dashed line). Based on our argu- 
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ments of Sec. |TT1 the charge spill-out is associated with 
a spatially decreasing D{z) [Fig. [TW b)]. This, in turn, 
modifies the interface potential barrier by producing a 
strong upward shift in energy of the conduction band 
edge from what one would have if D(z) were uniform and 
equal to the "physical" value Di. This implies that the 
charge spill out generally reduces the depolarizing field 
[the "pre-spill" estimate is sketched as a thick dashed line 
in Fig. HQf a)]. and hence overstabilizes the ferroelectric 
state. This is what one intuitively expects - population 
of the conduction band constitutes an additional channel 
for screening the polarization charge, and this cooperates 
with the metallic carriers of the electrode. This, however, 
contrasts with the conclusions of Ref. [H, where it was 
argued that charge leakage suppresses P by producing 
a ferroelectrically "dead" layer. These conclusions are 
based on the assumption that the physically measurable 
P is the average polarization, (P) , taken over the whole 
volume of the film. As the polarization is locally reduced 
near a pathological interface, charge spill-out indeed re- 
sults in a reduced (P) . 

Is it justified, though, to assume that (P) is the phys- 
ically relevant quantity in the capacitor? Does (P), in 
other words, refiect what is experimentally measured? 
In an experiment one measures the time integral of the 
transient current density, Aj, that flows through the ca- 
pacitor during the switching process. Aj does not relate 
to (P). Under the hypothesis that at least a portion of 
the film remains insulating throughout switching, it rig- 
orously follows from the modern theory of polarization^ 
that Aj = AD = 2\D\: D is the value of the (locally uni- 
form) electric displacement deep in the insulating region, 
indicated as D2 in Fig. [121 (We assume for simplicity 
that D = in the paraelectric reference state.) There- 
fore, observing that (P) is reduced upon charge leakage 
does not reflect the true physical effect of the pathologi- 
cal band alignment, which is an artificial enhancement of 
the spontaneous P via the reduction of the depolarizing 
field illustrated above. 

A large number of works^^"^"' have investigated the 
stability of PbTiOs-based capacitors, and it is impos- 
sible here to discuss in detail whether and how the above 
band-alignment issues might have affected each of them 
(for instance, regarding the polarization enhancements 
reported in Ref. l6l[) . We limit ourselves to observe that, 
due to the large spontaneous polarization of PbTiOa, the 
possible consequences of having a pathological ferroelec- 
tric state need to be taken seriously into account in the 
analysis, as we showed for the example of SrRuOa elec- 
trodes in Sec. El 



C. Transport properties in the tunneling regime 

Ferroelectric capacitors have been explored as poten- 
tial tunneling electroresistance devices;^ and many re- 
cent calculations focused on the calculation of the con- 
ductance by means of first-principles methods. Metallic- 



ity and spill-out of electrons is a serious potential issue 
in this context, as the calculated conductance can poten- 
tially be affected by the presence of space charge in the 
system, in a way which is difficult to predict. The recent 
work of Velev et aL^^ appears to be concerned by these 
issues, as it focuses on Ti02-terminated Pt/BaTiOs/Pt 
capacitors. Indeed, we find (see Sec. IIVB 2[) that this 
interface is problematic already in the centrosymmetric 
paraelectric case. While we have not explored the fer- 
roelectric regime in this system, based on the imperfect 
screening arguments of Sec. HIl (the lineup depends lin- 
early on P around the paraelectric reference phase) we 
expect the spill-out effect to become worse at least at one 
of the two interfaces when the capacitor is polarized. 

In fact, the metallicity of the ferroelectric film seems 
to be confirmed by the data presented by the authors: 
In Fig. 2(a-b) of Ref. [s^ the conduction band minimum 
(CBM) of the central BaTiOa cell appears to be degen- 
erate or lower than the Fermi level, and in Fig. 1 of the 
same paper the atomic displacements of the ferroelec- 
tric phase seems to be strongly asymmetric, consistent 
with our speculations. While we cannot draw a defini- 
tive conclusion (our computational setup slightly differs 
from that of Ref. [stI ). our analysis highlights the crucial 
importance of the band alignment issue, and the necessity 
of performing an adequate and convincing assessment of 
its impact on the results (e.g. the conductance) in each 
case. 



D. Interface magnetoelectric effects 

Magnetoelectricity is one of the emerging topics in ox- 
ide research. Despite the intense efforts, one of the main 
limiting factors still persists: bulk materials displaying a 
robust magnetoelectric effect are notoriously difficult to 
find. To work around this problem, several researchers 
have been looking for alternative solutions by exploring 
heterostructures and composite materials. An interface, 
due to its lower symmetry, might allow for physical re- 
sponse properties that are absent in the parent com- 
pounds. A promising route to interfacial ME coupling 
that has been proposed recentlj*^ is mediated by charge. 
The polarization of the ferroelectric (or dielectric) lattice 
produces a bound charge at the interface, that is screened 
by the carriers of the metal. If these carriers are spin- 
polarized, as in a ferromagnet, there will be a net change 
in the magnetization. 

It is easy to see that the band-alignment issues that we 
discuss in this work have direct and important implica- 
tions for the calculation of the carrier-mediated interface 
ME coefficient. In the pathological regime, the calcu- 
lated (magnetic) response will most likely be suppressed, 
as the spill-out charge, rather than the spin-polarized 
carriers in the electrode, will screen the applied bias po- 
tential (or the ferroelectric polarization). This specu- 
lation is directly relevant for interpreting the results of 
Yamauchi et al.^^ on BaTiOs films sandwiched between 
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Co2MnSi (Heusler alloy) electrodes. Depending on the 
termination, two qualitatively different behaviors were 
reported: the MnSi/Ti02 interface results in a patho- 
logical band alignment and a strongly non-homogeneous 
local polarization profile; conversely, neither is present in 
the capacitor with the other type of termination, which 
has symmetric Co/Ti02 interfaces. A very small magne- 
toelectric response was reported for the MnSi/Ti02 case 
(contrary to the Co/Ti02 case), in qualitative agreement 
with our arguments above. 

Other recent studiesj^i^ focusing on ME effects in 
thin Fe film deposited on ATiOs (A=Ba,Pb,Sr), also re- 
ported strongly non-uniform polarization profiles in the 
ferroelectric film (e.g. Fig. 3 of Ref . l68h . This suggests 
that also the ATiOa/Fe interface might be concerned by 
the band-alignment issues discussed in this work, with 
potential impact on the physical observables. Our anal- 
ysis tools should help clarify these issues in the above 
systems and in the Fe/BaTiOs/Fe capacitors of Ref. [gqI 
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FIG. 20: LDOS of selected Ti02 layers in the electrostatically 
doped LaAlOs/SrTiOa system. The insets are a blow up of 
the regions indicated by the small squares. The energy scale 
of the insets is comprised between -0.25 eV and 0.25 eV. 



F. Relationship to LaAlOa/SrTiOa 



E. Schottky barriers 

Direct calculations of Schottky barriers at 
metal/ferroelectric interfaces are, among the many 
useful physical properties of these junctions, those that 
are most directly affected by the issues we discuss here. 
The consequence of a pathological band alignment is 
that the estimated Schottky barrier is no longer a physi- 
cally meaningful interface property, but is infiuenced by 
macroscopic space-charge phenomena. 

A rather comprehensive work on the SrTiOs /transition 
metal interface was recently reported in Ref. [TO- Without 
going into too detailed an analysis of the results, we limit 
ourselves to noting that many of the reported p-type SBH 
for Ti02- or SrO-terminated interfaces are very close to, 
or sometimes well in excess of 1.8 eV. Considering that 
the LDA/GGA fundamental gap of SrTiOa is around 1.8 
eV, the actual n-type SBH of the calculation (i.e. not the 
value corrected with the experimental band gap) is close 
to zero or negative. Therefore, charge spill out is a con- 
crete and likely possibility for many of the investigated 
structures. 

Note that, contrary to the case of oxide electrodes, 
ideal interfaces between SrTiOs and simple metals tend 
to have a smaller Acff.--^ This implies that the effects 
of the electrostatic re-equilibration described in Sec. [ll] 
might be somewhat less dramatic, and the values of the 
self-consistent 0„ closer to (p^. This suggests that the 
trends and the conclusions reported in Ref. [t^ are likely 
to be robust with respect to the issues described in this 
work. However, a more detailed analysis would be cer- 
tainly interesting in order to assess their impact at the 
quantitative level. 



Many of the analysis tools developed in this work are 
not limited to ferroelectric capacitors, but can be readily 
extended to other systems where free-charge doping of a 
band insulator plays a central role. An excellent exam- 
ple, where the interpretation of the observed effects is still 
widely debated, is two dimensional conducting electron 
gas (2DEG) that forms at the polar LaAlOs/SrTiOa in- 
terface 1^ A central problem is the determimation of the 
physical effects governing the confinement and equilib- 
rium distribution of the 2DEG. Some authors^^ propose 
a mechanism for the confinement of the gas based on 
the formation of metal induced gap states (MIGS) in the 
band gap of SrTiOa. Other authors^ however, explain 
the experimental observations in terms of a semiclassical 
Thomas-Fermi model that is analogous to that described 
in Sec. HIl and where the MIGS are completely absent. 
Answering the question of whether the MIGS play an im- 
portant role in this system involves a careful analysis of 
the local electronic properties, and more specifically of 
the LDOSiiS In this sense, the methodology discussed in 
Sec, nil A~3l appears ideally suited to clarifying this issue. 

We base our analysis on the calculations done in 
Ref. iM, with a 24-cell thick SrTiOa slab and a 3-cell 
LaAlOs overlayer. (This calculation was performed with 
a 12 X 12 Monkhorst-Pack sampling of the surface Bril- 
louin zone, and with a Gaussian smearing of 0.1 eV; full 
details of the computational parameters are reported in 
Ref. Q.) The boundary conditions are set to I?sto=0, 
DhAO = "6/25*, and are equivalent to those of the sym- 
metric superlattice used by Janicka et al^ In Fig. [20lwe 
show the LDOS corresponding to the Ti02 layers number 
15 (curve a), 10 (b) and 5 (c), where layer 1 is adjacent to 
the LAO interface. On top of each curve we superimpose 
the bulk Ti02 LDOS, that we align with the supercell 
LDOS by matching the semicore Ti(3s) peak at ^ —57.5 
eV. (As in earlier works the 0(2s) peak was used as a 
reference, we also show the 0(2s)-derived feature, which 
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is located at about -18 eV.) The matching is excellent 
in all cases, especially in the layers lying furthest from 
the interface where the effect of the structural distor- 
tions and free charge are less pronounced. (Note that 
we performed the bulk calculation with a /c-point mesh 
that accurately matches the one used in the supercell cal- 
culation. Also, in the construction of the LDOS curves 
we used the same Gaussian smearing function of width 
0.1 eV, corresponding to the smearing used to relax the 
self-consistent ground state of the supercell structure.) 
In the insets we show a blow-up of the conduction band 
edge, which goes below the Fermi energy in agreement 
with the semiclassical arguments of Ref. [tI and of our 
Sec. HH Clearly, our plots do not show any evidence for 
MIGS in the energy gap, contrary to the conclusions of 
Janicka et alJ"^ 

To reconcile this discrepancy, we can speculate that the 
LDOS curves presented in Ref. might have been con- 
structed with a substantially larger smearing width than 
ours, and this might have precluded an accurate iden- 
tification of the band edges. We believe that the tech- 
nique presented here (of superimposing an appropriately 
constructed bulk LDOS on top of the supercell curves) 
provides a very practical means of minimizing systematic 
errors in the analysis of the results. 



VII. CONCLUSIONS 

Due to its accuracy and efficiency, density functional 
theory has emerged as the method of choice for study- 
ing ferroelectric oxides from first-principles. This pre- 
dominance has been reinforced since the early 1990s by 
the many successes achieved in the determination of the 
structural, energetic, piezoelectric, and dielectric proper- 
ties at the bulk level. In the last few years, those efforts 
have evolved to address the behaviour of the functional 
properties in thin films and superlattices, including in 
same cases (for instance, in the study of ferroelectric ca- 
pacitors) the presence of metal/insulator interfaces. 

For a reliable prediction of the functional properties 
of these devices, the atomic displacements, distortions of 
the unit cell, the electronic structure and the band gap 
have to be accurately described simultaneously. However, 
the proper DFT treatment of such interfaces is compli- 
cated by the so-called "band-gap problem" , which might 
produce a pathological alignment between the Fermi level 
of the metal and the conduction band of the insulator, 
thus precluding explicit DFT investigation of many sys- 
tems of practical interest. In this work we provide useful 
guidelines to identify such a pathological scenario in a 
calculation by examining its main physical consequences: 
(i) an inhomogeneous polar distortion propagating into 
the bulk of the film, (ii) the film becoming partially or to- 
tally metallic due to a non- vanishing free charge, and (iii) 
the local conduction band edge crossing the Fermi level. 
The above three effects are intimately linked, and should 
be considered as potential artifacts of the aforementioned 



band-gap problem. Whenever one of these "alarm fiags" 
is raised in a calculation, the results should be examined 
with great caution. 

A route to overcoming this limitation involves correct- 
ing the LDA/GGA bandgap while preserving the excel- 
lent accuracy of these functionals in the prediction of 
ground-state properties. Traditional methods to increase 
the band gap of insulators, like the inclusion of a Hubbard 
U term in the Hamiltonian, are not satisfactory in the 
case of a ferroelectric capacitor with B-cation driven fer- 
roelectricity. A more promising avenue has been recently 
opened by Bile et alJ^ and Wahl and coworkers, using 
the so-called "hybrid" functionals that combine Hartree- 
Fock exchange and DFT. In particular the Bl-WC func- 
tional proposed in Ref. 't^ has been shown to provide 
good structural, electronic and ferroelectric properties as 
compared to experimental data for BaTiOs and PbTiOa; 
verifying the accuracy of Bl-WC in interface studies will 
be an interesting subject for future research. Unfortu- 
nately, the price to pay for this accuracy is the substan- 
tially higher computational cost of Bl-WC compared to 
LDA/GGA. 

In addition to the purely technical issues, our work also 
opens interesting avenues regarding fundamental physical 
concepts. For example, ferroelectricity is usually under- 
stood within the modern theory of polarization, which 
is only applicable in the absence of conduction electrons 
(i.e. in pure insulators at zero electronic temperature). 
It is an important fundamental question, therefore, to as- 
sess whether our understanding of ferroelectrics in terms 
of bound charges, polarization and macroscopic electrical 
quantities still applies (and to what extent) in a regime 
where a sizable amount of space charge is present in the 
system. This issue is of crucial importance also for other 
systems, e.g. electrostatically doped perovskites, which 
bear many analogies to the physical mechanisms dis- 
cussed in this work. The first-principles-based modeling 
approach proposed in Ref. [sj appears to be a promising 
route to further exploring this interesting topic. 
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etot 


too 


Xtot/xion 


BaTiOs 


-48.87 


6.48 


0.90 


PbTiOs 


-96.54 


8.33 


0.93 


KNbOa 


-34.92 


6.27 


0.87 



TABLE III: Values of the susceptibilities x and scaling factors 
Xtot/xion for the ferroelectric materials considered in this 
work. 



of Cantabria, and at CESGA. 



Appendix A: Local polarization via Born effective 
charges 

In this Appendix we discuss the approach, used in sev- 
eral parts in this manuscript and ubiquitously in the re- 
cent literature, of associating the local value of the "effec- 
tive" polarization (i.e. the induced P with respect to the 
reference centrosymmetric configuration^) in capacitor 
heterostructures with an approximate formula, based on 
the Born effective charges, Z* . In particular, we provide 
formal justification for an improved formula, still based 
on the Z* , that we introduced in this work, and we al- 
ready mentioned in Sec. IIIIB 31 

Recall the definition of the approximate effective po- 
larization in terms of the Born effective charges in a bulk 
sohd. 



— ^ ^ Z^Raz 
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■are Z 
Berry phase 
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Reduced electric displacement d 

FIG. 21: Polarization P in a BaTiOa film and PbTiOa bulk as 
a function of the reduced electric displacement field d = DS. 
Data are taken from Ref. [20l (see Section III.C.l) and Ref . l3^. 
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where ctot is the total dielectric constant of the insulator 
(relative to the vacuum permittivity eo). Substituting 
Eq. (jA3| and Eq. (jAU into Eq. (|A2l) 



(A5) 



XION 
eoCTOT 



The same kind of arguments applied to the total po- 
larization yield 



It is easy to verify that the layer-resolved expression 
of Eq. ([5^ reduces to P^ in the case of a periodic crys- 
tal, where P^ is a constant function of the layer index j. 
P^ does not reduce to the "correct" polarization P{D) 
at any value of D, as it does not take into account the 
additional polarization of the electronic cloud due to the 
internal field £{D) (recall that the Born effective charges 
are defined under the condition of zero macroscopic elec- 
tric fieldJ^) 

Taking the Taylor expansion of the polarization as a 
function of D (we assume for simplicity that D, P and P^ 
all vanish in the reference centrosymmetric structure) , we 
can write 



P^{D) ^ 



dP^ 

Id 



-D 



dP^ dE 
d£ dD' 



-D 



(A2) 



For small values of I?, we can truncate the previous ex- 
pansion at the linear order term. Now, by definition 



dP^ 
d£ 



XlON, 



(A3) 



where xion is the lattice-mediated susceptibility, and 



P{D) - D 



Xtot 
coEtot 



(A6) 



where xtot is the sum of the lattice-mediated suscepti- 
bility, xiON, and the purely electronic (frozen- ion) sus- 
ceptibility, Xoo- Note that xion is not bound to be pos- 
itive. In a ferroelectric material, for example, the cen- 
trosymmetric reference structure is unstable and there- 
fore yields a negative Xion (and hence ctot); as dis- 
cussed in Ref. |32| . The present derivation is general and 
encompasses those cases. 

From the above considerations it immediately follows 
that an estimate of the total polarization, which is exact 
in the linear limit, can be given as 



P{D) 



Xion 



-P^{D). 



(A7) 



This is essentially Eq. ((33)) . In practice, xion and Xoo 
are calculated in the reference phase according to the 
standard definitions,— 



Xion = eo(eTOT — foe) — 



E 



(A8) 
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where Mq is a unit mass, are the normal mode charges 
and are the eigenvalues of the dynamical matrix, and 



Xoo = eo(eoo - 1), 



dD 



fixed — ions 



(A9) 



The values of these physical constants that are relevant 
for the results presented in this manuscript are reported 
in Table Uni 

We proceed in the following to test this approxima- 
tion on two representative bulk ferroelectric materials, 
PbTiOa and BaTiOa. We take the relevant data (lin- 
ear susceptibilities. Born charges and relaxed structures 
as a function of D) from the calculations of Ref. HH and 
Ref.H^. Note that the BaTiOs calculation was performed 
at a fixed value of the in-plane lattice parameter (indi- 
cated as "film" in the figure) while in the PbTiOa cal- 
culation both a and c parameters were relaxed for each 
value of D. The results are presented in Fig. [5TJ In 
both cases, the "bare" value is systematically overes- 
timated compared to the Berry-phase polarization. With 
the correction described above, i.e. by rescaling all val- 
ues by the factor xtot/xion, the approximate value of 
P accurately matches the Berry-phase one. The accu- 
racy is surprisingly good in BaTiOs, where the maximum 
deviation is of the order of 1%. In PbTiOa, for large val- 
ues of d the rescaled-Z* value of P presents significant 
deviations. Note that these deviations mostly concern 
values of d that are larger than that of the ferroelectric 
ground state {d ~ 0.74), and therefore are not of concern 
in this manuscript. We ascribe these deviations to the 
field-induced structural transition that was described in 
Ref. m. 

In conclusion, this simple rescaling factor appears to 
be an effective way to obtain a relatively accurate value 
of the local P in heterostructure calculations, based only 
on the local atomic positions and a few ingredients that 
can be easily computed in the bulk reference structure. 
From the results of our tests, we expect the agreement to 
be best in cases where the polarization is small (closer to 
the linear limit where the approximation becomes exact). 
Furthermore, cases where the ferroelectric polarization 
can be represented in terms of a single "soft mode" such 
as BaTiOa seem to work better than cases, like PbTiOa, 
where significant mode mixing and non-trivial structural 
transitions occur at higher D values. 



Appendix B: Convolutions and energy smearing of 
the local density of states 

1. Convolutions 

Convolution is a mathematical operation on two func- 
tions / and g, producing a third function that is typically 
viewed as a modified version of one of the original func- 
tions. For the purpose of the present notes, it is useful to 



think of / as a data curve containing the relevant physi- 
cal information, and 5' as a rapidly-decaying "smoothing" 
function that produces a local weighted average of /. We 
define the convolution of / and g, f * g, as the following 
integral transform. 



-\-oo 



f{v)g{x-v)dy 



(Bl) 



Convolutions have many properties, including commu- 
tativity and associativity. Furthermore, the Dirac delta 
can be thought as the identity under the convoluton op- 
eration. 



(B2) 



and under certain assumptions an inverse operation can 
also be defined. In other words, the set of invertible dis- 
tributions forms an abelian group under the convolution. 

A particularly useful property holds in relationship to 
the Fourier transform. 



T{f*g) = k-F{.f)-F{g) 



(B3) 



where J^{f) denotes the Fourier transform of /, and k is 
a constant that depends on the normalization convention 
for the Fourier transform. Thus, in reciprocal space the 
convolution becomes a simple product. This naturally 
provides an efficient convolution algorithm; the workload 
is reduced from ©(iV^) to 0[N\og{N)]. 



2. Local density of states 

In this work we use [Eq. ([H])] the following formula to 
compute the smeared local density of states (LDOS), 

p(r, E)^Y.^^ IV'«k(r)|' g{E ~ E^k). (B4) 

nk 

We shall see that this is indeed a convolution. We first get 
rid of the spatial cordinates. To this end, it is customary 
to integrate the LDOS in real space over a given volume 
V, 

Pv{E) = Y,WkPnk{V)g{E - S„k), (B5) 



where 



P„k(V^) 



rfV|V'„k(r)r 



(B6) 



(Note that sometimes it might be more convenient to use 
a projected density of states, rather than a local density 
of states. In such cases it is sufficient to replace the real- 
space integral in the above equation with an appropriate 
sum over angular momentum components. The follow- 
ing discussion remains unchanged.) Now the LDOS is a 
function of a single energy variable. If we write 



fv{E) - ^ii;kP„k(X^)<5(S- S„k), 



(B7) 



nk 
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we can easily see that pv = fv * g- This leads to a 
simple reciprocal-space expression. We first define an en- 
ergy window, [i?iow, -Bhigh], that contains the entire eigen- 
value spectrum Enk. We actually take a window which 
is slightly larger, where this "slightly" depends on the 
decay properties of g, 

Blow = min(£'„k) - e, £^high = max(£;„k) + £■ (B8) 

The width of this window is -Ehigh — -Eiow = Ai?. We 
represent pv{E) in reciprocal space as a discrete Fourier 
transform, 



(B9) 
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where u — 21m/ /S.E and n is an integer. By using 
Eq. ([B3| we have 



= l\E ■ fv{uj) ■ g{uj). 



(BIO) 



The Fourier transform of a Dirac delta centered in the 
origin is a constant. Eq. (jBlOp then decomposes the local 
density of states into a structure factor^ 



1 



/yM-^E^kP„k(l^)e-*'^"-, (Bll) 

Tik 

and a form factor g{uj). Obviously, this formulation is 
only convenient if the function g has a fast decay in both 
real and reciprocal space, so that the sum in Eq. (jB9|) 
can be truncated. This is indeed the case for the most 
widely used smoothing functions g, as we shall see in the 
following. 



3. Gaussian vs. Fermi-Dirac smearing 

The Gaussian smearing (G) and the Fermi-Dirac (FD) 
smearing are by far the most popular choices for the occu- 
pation function in first-principles calculations of metallic 
systems. If we define the occupation function / as the 
integral of a "kernel" function g, 



fiE) = 1 



g{x)dx, 



(B12) 



one can verify that the Gaussian or Fermi-Dirac occupa- 
tion are, respectively, reproduced by the following choices 
of 9, 



gcix) 



gpuix) = 



/ira 



2 + e=^/'^ -I- e-="/'^' 



(B13) 
(B14) 



where a is the smearing energy [these correspond to 
Eq. (|23a|) and Eq. ()23bp ]. It is easy to see that, by com- 
bining Eq. (|B13|) or Eq. (|B14p with Eq. (|B12p one obtains 
the standard definitions of the occupation function (we 



FIG. 22: (a) Gaussian (cr = 0.15 eV) and Fermi-Dirac (a = 
0.075 eV) occupation functions, (b) Kernel of the occupation 
functions as defined in tlie text, (c-d) Fourier transform of the 
smearing kernels g, assuming an energy window of [—1, 1]. 



assume that the complementary error function, erfc, val- 
ues 2 at — oo), 



foix) = 5erfc(a;/CT), 
/FD(a;) = gx/Hr^i - 



(B15) 
(B16) 



It is useful to spell out the explicit formulas for the 
Fourier transforms of both smearing functions. 



5g(w) 

gpui^) 



AE ' 
AE sinh(7rw(T) 



(B17) 
(B18) 



Note that the above formulas are normalized according to 
the conventions on the Fourier transforms that we used 
in the previous section. The functions / and g defined 
above are shown in Fig. [511 Note that a different choice 
of a was used in the Fermi-Dirac and in the Gaussian 
case. A FD distribution is roughly equivalent to a G 
distribution with a a value that is twice as large. 

In the main text and here we have assumed that it is a 
good idea to use the same g kernel in the calculation and 
in the construction of the LDOS. We shall substantiate 
this point in the following Section. 



4. On the optimal choice of g 

In many cases, the specific choice of the g function 
to be used in Eq. (IB4p is largely arbitrary. Typically, 
the goal is to filter out the unphysical wiggles due to 
the discretization of the fc-mesh, but at the same time 
to preserve the main physical features, without blurring 
them out completely. This calls for a smearing function 
that is neither too sharp nor too broad. Since a "slightly 
too broad" or a "slightly too sharp" smearing function 
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Energy (eV) Energy (eV) 

FIG. 23: Left: Fermi-Dirac occupation function, identical to 
that of Fig.[22ja) (solid curve); hypothetical orbital located at 
an energy of 0.15 eV above the Fermi level (dashed line); the 
thermal occupation of this state yields a total charge of 0.119 
electrons (red dot). Right: density of states corresponding to 
the single isolated orbital at an energy of 0.15 eV above the 
Fermi level, smeared by using the qfd kernel of Fig. [22jb); 
the integral of the DOS up to the Fermi level (shaded area) 
yields the exact same charge of 0.119 electrons. 

usually does not influence the physical conclusions, in 
many cases one has the freedom of choosing whatever 
yields the clearest visual aid to support the discussion. 

There are cases, however, where this choice is not just 
a matter of aesthetics, and using the "wrong" g function 
can qualitatively and quantitatively influence the inter- 
pretation of the results. More specifically, the issue con- 
cerns cases where the analysis of the LDOS (or DOS or 
PDOS) is used to detect and quantify the population of 
orbitals that lie close in energy to the Fermi level. As 
we focus on charge spill-out phenomena that concern the 
conduction band of a dielectric/ferroelectric film in con- 
tact with a metallic electrode, this is a central point of our 
work. The problem is most easily appreciated by looking 
at the left panel of Fig. [221 There is a single orbital lying 
at an energy of 0.15 eV above the Fermi level. As this or- 
bital lies above the Fermi level, one might be tempted to 



think that the orbital is empty, and that charge spill-out 
does not occur at all. However, calculations in metallic 
systems are routinely performed by using an occupation 
function that is artificially broadened, in order to improve 
convergence of the ground-state properties; in Fig. 1231 we 
assume a Fermi-Dirac occupation with a fictitious elec- 
tronic temperature of 0.075 eV. It is easy to see that with 
such an occupation function, the orbital lying at 0.15 eV 
won't be empty, but will be "thermally" populated by 
tail of the Fermi-Dirac distribution. The final result is a 
charge transfer of 0.119 electrons into this orbital. 

Now, is there a "right" way to construct the DOS 
curve, such that the above-mentioned charge transfer 
could be qualitatively and quantitatively inferred from 
the DOS, without knowing any further detail of the cal- 
culation? The answer is yes, and consists in constructing 
the DOS by using as broadening g function which is con- 
sistent with the occupation function used by the code. 
In this case, this is gpD, with a a identical to that used 
to calculate the electronic ground state. To demonstrate 
this point, we plot in the right panel of Fig.[221the DOS of 
this isolated orbital at 0.15 eV, appropriately convoluted 
with (7fd- Eq. (jB12[) guarantees that, by doing this, one 
recovers the very intuitive result that the total amount of 
electron charge, Q, present in the volume V (over which 
the LDOS was integrated) exactly corresponds to the in- 
tegral of the DOS up to the Fermi level, 

pv{E)dE. (B19) 

-oo 

Then, a simple look at the DOS curve is sufficient to 
ascertain whether a significant transfer of charge has oc- 
curred into a specific group of bands. As this rigorous 
sum rule can be very practical in the analysis of the re- 
sults, we encourage a systematic use of the "internally 
consistent" LDOS construction described above. 



1 J. F. Scott and C. A. P. de Araujo, Science 246, 1400 
(1989). 

^ O. Auciello, J. F. Scott, and R. Ramesh, Physics Today 

July, 22 (1998). 
^ J. F. Scott, Ferroelectric memories (Springer- Verlag, 

2000). 

* M. Dawber, K. M. Rabe, and J. F. Scott, Rev. Mod. Phys. 
77, 1083 (2005). 

^ J. F. Scott, J. Phys.: Condens. Matter 18, R361 (2006). 
^ J. F. Scott, Science 315, 954 (2007). 
M. Fiebig, J. Phys. D: Appl. Phys. 38, R123 (2005). 

* W. Eerenstein, N. D. Mathur, and J. F. Scott, Nature 
(London) 442, 759 (2006). 

^ C. H. Ahn, J. Mannhart, and J.-M. Triscone, Nature (Lon- 
don) 424, 1015 (2003). 

E. Y. Tsymbal and H. Kohlstedt, Science 313, 181 (2006). 
P. Maksymovych, S. Jesse, P. Yu, R. Ramesh, A. Baddorf, 
and S. V. Kalinin, Science 324, 1421 (2009). 
V. Garcia, S. Fusil, K. Bouzehouane, S. Enouz-Vedrenne, 



N. Mathur, A. Barthelemy, and M. Bibes, Nature (London) 
460, 81 (2009). 

V. Garcia, M. Bibes, L. Bocher, S. Valencia, F. Kronast, 
A. Crassous, X. Moya, S. Enouz-Vedrenne, A. Gloter, 
D. Imhoff, et al.. Science 327, 1106 (2010). 
P. Ghosez and J. Junquera, in Handbook of Theoretical and 
Computational Nanotechnology, edited by M. Rieth and 
W. Schommers (American Scientific Publisher, Stevenson 
Ranch, CA, 2006). 

J. Junquera and P. Ghosez, J. Comput. Theor. Nanosci. 5, 
2071 (2008). 

I. Souza, J. Ifiiguez, and D. Vanderbilt, Phys. Rev. Lett. 
89, 117602 (2002). 

P. Umari and A. Pasquarello, Phys. Rev. Lett. 89, 157602 
(2002). 

I. Souza, J. Ifiiguez, and D. Vanderbilt, Phys. Rev. B 69, 
085106 (2004). 

M. Stengel and N. A. Spaldin, Phys. Rev. B 75, 205121 



29 



(2007). 

^° M. Stengel, D. Vanderbilt, and N. A. Spaldin, Phys. Rev 
B 80, 224110 (2009). 

J. P. Perdew and M. Levy, Phys. Rev. Lett. 51, 1884 
(1983). 

L. J. Sham and M. Schliiter, Phys. Rev. Lett. 51, 1888 
(1983). 

M. Stengel and N. A. Spaldin, Nature (London) 443, 679 
(2006). 

J. Junquera and P. Ghosez, Nature (London) 422, 506 
(2003). 

M. Y. Zhuravlev, R. F. Sabirianov, S. S. Jaswal, and E. Y. 
Tsymbal, Phys. Rev. Lett. 94, 246802 (2005). 
V. Heine, Phys. Rev. 138, A1689 (1965). 
S. Okamoto and A. J. Millis, Nature 428, 630 (2004). 
^® A. Baldereschi, S. Baroni, and R. Resta, Phys. Rev. Lett. 
61, 734 (1988). 

L. Colombo, R. Resta, and S. Baroni, Phys. Rev. B 44, 
5572 (1991). 

^° J. Junquera, M. H. Cohen, and K. M. Rabe, J. Phys.: 
Condens. Matter 19, 213203 (2007). 

A. Franciosi and C. G. V. de Walle, Surface Science Re- 
ports 25, 1 (1996). 

M. Stengel, N. A. Spaldin, and D. Vanderbilt, Nature 
Physics 5, 304 (2009). 

M. Stengel, D. Vanderbilt, and N. A. Spaldin, Nature Ma- 
terials 8, 392 (2009). 
M. Stengel, arXiv: 1012.2468 (2010). 

M. Peressi, N. Binggeh, and A. Baldereschi, J. Phys. D: 
Appl. Phys. 31, 1273 (1998). 
^® In practical simulations, the origin of the fc-point grid may 
be displaced from fc = in order to decrease the number of 
inequivalent fc-points.— This shift usually prevents the 
appearance of high-symmetry points from the list of k- 
points used during the self-consistent procedure or in the 
calculations of the density of states. 

The band gap of both BaTiOs and PbTiOs is indirect, 
with the top of the valence band located at R in BaTiOa 
and at X in PbTiOa, and the bottom of the conduction 
band at T in both materials. 

^® K. T. Delaney, N. A. Spaldin, and C. G. Van de Walle, 
Phys. Rev. B 81, 165312 (2010). 

3^ J. Bardeen and W. Shockley, Phys. Rev. 80, 72 (1950). 
C. J. Fall, N. Binggeh, and A. Baldereschi, J. Phys.: Con- 
dens. Matter 11, 2689 (1999). 

A. Grigoriev, R. Sichel, H. N. Lee, E. C. Landahl, 

B. Adams, E. M. Dufresne, and P. G. Evans, Phys. Rev. 
Lett. 100, 027604 (2008). 

X. Wu, O. Dieguez, K. M. Rabe, and D. Vanderbilt, Phys. 
Rev. Lett. 97, 107602 (2006). 
■'^ E. D. Murray and D. Vanderbilt, Phys. Rev. B 79, 100102 
(2009). 

N. Marzari and D. Vanderbilt, Phys. Rev. B 56, 12847 
(1997). 

J. B. Neaton and K. M. Rabe, Appl. Phys. Lett. 82, 1586 
(2003). 

■'^ M. Stengel and D. Vanderbilt, Phys. Rev B 80, 241103 
(2009). 

'^'^ J. M. Soler, E. Artacho, J. D. Gale, A. Garcia, J. Junquera, 
P. Ordejon, and D. Sanchez- Portal, J. Phys.: Condens. 
Matter 14, 2745 (2002). 

*® M. Payne, M. Teter, D. Allan, T. Arias, and J. Joannopou- 



los. Rev. Mod. Phys. 64, 1045 (1992). 
"'^ P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 

H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 
(1976). 

" J. Moreno and J. M. Soler, Phys. Rev. B 45, 13891 (1992). 
L. Kleinman and D. M. Bylander, Phys. Rev. Lett. 48, 
1425 (1982). 

N. Troullier and J. L. Martins, Phys. Rev. B. 43, 1993 
(1991). 

J. Junquera, M. Zimmer, P. Ordejon, and P. Ghosez, Phys. 
Rev. B 67, 155327 (2003). 
^5 C. H. Peng, J. F. Chang, and S. Desu, MRS Symp. Proc. 
243, 12 (1992). 

C.-G. Duan, R. F. Sabirianov, W.-N. Mei, S. S. Jaswal, 

and E. Y. Tsymbal, Nano Letters 6, 483 (2006). 
" J. P. Velev, C.-G. Duan, K. D. Belashchenko, S. S. Jaswal, 

and E. Y. Tsymbal, Phys. Rev. Lett. 98, 137201 (2007). 

Y. Wang, M. K. Niranjan, K. Janicka, J. P. Velev, M. Y. 

Zhuravlev, S. S. Jaswal, and E. Y. Tsymbal, Phys. Rev. B 

82, 094114 (2010). 
5^ J. D. Burton and E. Y. Tsymbal, Phys. Rev. B 82, 161407 

(2010). 

R. Resta and D. Vanderbilt, in Physics of Ferroelectrics: A 
Modem Perspective, edited by K. M. Rabe, C. H. Ahn, and 
J.-M. Triscone (Springer- Verlag, Berhn Heidelberg, 2007). 
^'^ N. Sai, A. M. Kolpak, and A. M. Rappe, Phys. Rev. B 72, 
020101 (R) (2005). 

Y. Umeno, B. Meyer, C. Elsasser, and P. Gumbsch, Phys. 
Rev. B 74, 060101 (R) (2006). 

W. A. Al-Saidi and A. M. Rappe, Phys. Rev. B 82, 155304 
(2010). 

^'^ Y. Umeno, J. M. Albina, B. Meyer, and C. Elsasser, Phys. 
Rev. B 80, 205122 (2009). 

J. M. Rondinelli, M. Stengel, and N. A. Spaldin, Nature 
Nanotechnology 3, 46 (2008). 

K. Yamauchi, B. Sanyal, and S. Picozzi, Appl. Phys. Lett. 
91, 062506 (2007). 
"'^ M. Fechner, I. V. Maznichenko, S. Ostanin, A. Ernst, 
J. Henk, P. Bruno, and I. Mertig, Phys. Rev. B 78, 212406 
(2008). 

M. Fechner, I. V. Maznichenko, S. Ostanin, A. Ernst, 

J. Henk, and I. Mertig, Phys. Status Sohdi B 247, 

16001607 (2010). 
®^ C.-G. Duan, S. S. Jaswal, and E. Y. Tsymbal, Phys. Rev. 

Lett. 97, 047201 (2006). 
™ M. Mrovec, J.-M. Albina, B. Meyer, and C. Elsasser, Phys. 

Rev. B 79, 245121 (2009). 

A. Ohtomo and H. Y. Hwang, Nature 427, 423 (2004). 

K. Janicka, J. P. Velev, and E. Y. Tsymbal, Phys. Rev. 

Lett. 102, 106803 (2009). 
'^^ O. Copie, V. Garcia, C. Bodefeld, C. Carretero, M. Bibes, 

G. Herranz, E. Jacquet, J.-L. Maurice, B. Vinter, S. Fusil, 

et al., Phys. Rev. Lett. 102, 216804 (2009). 
'^'^ D. I. Bile, R. Orlando, R. Shaltaf, G.-M. Rignanese, 

J. Ifiiguez, and P. Ghosez, Phys. Rev. B 77, 165107 (2008). 
'^^ R. Wahl, D. Vogtenhuber, and G. Kresse, Phys. Rev. B 

78, 104116 (2008). 
'^^ P. Ghosez, J.-P. Michenaud, and X. Gonze, Phys. Rev. B 

58, 6224 (1998). 
■^■^ X. Gonze and C. Lee, Phys. Rev. B 55, 10355 (1997). 



